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ABSTRACT 


We ventured to develope a complete software package 

for fault tree analysis of any system be it mechanical, 

electrical or nuclear. In doing so attempt was made to 

improve upon the existing computer codes in terms of ease to 

use, CPU time, efficiency and software facilities, A 

particular attention was paid to enlarge the scope of the 

software. For example, almost all kinds of gates have been 

allowed using fault tree modification program. In fault 

tree quantification, ’M’ out of *N* type of system has been 

allowed and a new formulation has been dbne using Fussell's 

approximation. Save fault tree construction, all other 

processes in the fault tree analysis has been automated 

including the electrical systems associated in nuclear systems, 

W —7 

We got unavailability at the end of 40 years ^0.5x10 for 
SNR-300 primary loop which is satisfactoiry according r to 
10 CFR 100 criterions,* We also established that increase 
of redundancies improves the system but some componenta being 
more critical than others. Our computer codes can, hende, 
be of immense help in system design. 



CHAPTER 1 


INTRODUCTION 


The present work is devoted to generate a complete 
software package for reliability analysis of any system 
using fault tree method in a most general form. Also 
discussed is the application of the package to primary 
coolant loop with emergency core cooling of a LMFBR for 
which model chosen is that of SNR-300 (A Geiman 300 MW 
Fast Ereader Reactor) • This package can be compared 
favourably with earlier similar packages in terms of its 
facilities, methodology* efficiency and CPU time on a 
computer. 

The first program finds minimal cutsets for a 

OR 

fault tree containing upto lOO^gates or 100 basic events. 
This number could be increased easily by increasing dimen- 
Sion of the concerned variables, MOCUS, on the other hand, 
can be used to find minimal cuts for upto 20 gates in a 
given tree. In this new program PCOMCP all basic events 
are coded by prime numbers and all logical operations 
and minimization processes are carried out by simple 
arithmetic operations. This reduces the storage require- 
ment greatly since no character reading or logical 
operations are needed, A new method of providing input 
has been used in which number of lines or cards is no 
more than the number of basic events; thus simplifying 
the input process a great deal. In most of other programs 
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input IS given gate by gate .Thus, number of input cards 
equals number of gates and hence a lengthy process. As 
in MOCUS, the maximum order of cutsets can also be specified 
in the input data. 

Second program quantifies the fault tree, i.e, 
finds top event unavailability (probability) and other 
related parameters. This program is based on KITT G3Q 
formulation but there is an improvement. The present program 
takes care of M out of N systems for basic events and top 
event. It is much easier to specify the values of M & N 
rather than modifying fault tree itself to take into 
account such systems. This kind of systems are mostly 
encountered due to large redundancy in nuclear systems to 
improve reliability. As in KITT1 and KITT2 bracketing 
procedure has been used to find maximum and minimum top 
event unavailability. In present program, at some places 
FUSSEL's approximation has been used when such parameter are 
not accurately desirable and formulation for exact value is 
otherwise difficult. Data could be provided both in terms 
of Lambda and Meu or Q (unavailability) and w (unconditional 
failure intensity). In latter case by bisection method erf 
iteration , Lambda and Meu are obtained. There are three 
options for providing accuracy desired. Option 0 
take all cutsets into account. Option 1 =>take maximum 
order of cutsets ( specified). Option > 1 =£> specify the 
% accuracy desired. Integration uses Simpson's rule in 
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which deciional place accuracy is specified. It is obvious 
that program is most versatile and flexible in use. 

Failure data and typical repair data for basic 

components have been obtained from WASH - 1400 report[.V^3 

Nuclear Engg. & Design (1984 V 81-82, NED 1984 V 83), 

n't! 

Reliability technology by Green and Bourne^and other 

sources. From the analysis an unavailability of 
-7 

0,5x10 has been obtained at the end of 40 year and it is 
seen that with time it decreases but not at a large rate 
for the primary coolant loop of LMFBR, This low unavail- 
ability has been obtained in LMFBR by using sufficient 
redundancies and standby systems. 

Third program finds minimal cutsets for an 
electrical systems which contain large number of feed- 
backs, interconnection, and are represented by a ckt, 
graph rather than a fault tree. Entirely different approach 
IS needed in this case. Electrical systems are integral 
part of nuclear power plant hence the interest. The 
program finds basic minimal path and the combinations 
which break these paths constitute the minimal cutsets 
[Chapter 6 ] , 

Last program is a fault tree modification program. 

It modifies the fault tree to take care of XDR, NOT, NAND, 
NOR gates. It utilises the gate equivalence for XOR and 
De Morgan’s law for others. In some cases, level of the 
tree increases which makes the program slightly complex. 
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This kind of facility is not available in MOCUS or PREP- 
KITT, 'M* out of 'N* kinds of gates are taken care of 
at quantification level. Output of this program is a fault 
tree containing only AND and OR gates. Some of the basic 
events occur in negated form. It first determines the 
sub-branches which necessitates increase in the level of 
tree. All the basic event entries are modifies accordingly, 
GATE equivalence and De Morgan’s law suggest the rest of 
the procedure. 



CHAPTER 2 


SYSTEM DESCRIPT lOM AND MODELLING 


The system we chose in this study for reliability 
analysis was primary sodium loop and emergency core cooling 
system of SNR— 300 which is also called Kalkar nuclear power 
station (a joint German, Belgium and Dutch project). This 
design is very similar for many fast breeders and as we 
will see the main features are also valid for Clinch-River 
Breeder reactor plant and other fast breeder reactors. 

2.1 SNR 300 Details (Only Primary & Intermediate Loop)r } 

The three parallel loops of the heat transfer system 
(HTS) transfer the heat produced in the reactor via inter- 
mediate heat exchangers and steam generators to the 
electricity producing system. They are furthermore used for 
the decay heat removal with the pump running at 5% of 
their nominal speed. Even with two loops out of service t 
natural convection in the third loop can remove the decay 
heat without reaching boiling at the core outlet. Even 
in cases leading to a failure of all main loops, the decay 
heat removal can take place through six emergency heat 
exctengprs operating in parallel and located inside the 
reactor tank. 
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The pipes and components of heat transfer systems with the 
exception of the steam generator are using unstabilized 
austenitic stainless steel X8 Cr Ni 18t1 (W1,4948), 

Each primary loop contains one circulating pump, a bank 
of three intermediate heat exchangers, one sodium flow meter 
and one throttle valve to reduce the coolant flow following 
a reactor shut down. 

The coolant circulating pumps are positioned in 
the hot leg to provide a sufficient suction head at a low 
excess pressure prevailing in the upper part of the reactor 
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tank. The pumps are of variable speed, single stage, 
single-flow, radial type, arranged vertically with free 
sodium surface and argon cover gas atmosphere. 

Heat exchangers are of straight tube design with 
a floating head at the lower end mth the primary sodium 
flowing on the shell side, an arrangement that allows the 
tube bundle to be withdrawn without cutting the primary 
loop pipes. 

The pipe connecting the reactor vessel to the pump 
has a diameter 600 mm which reduces to 500 mm for all 
other main pipes, 

A butterfly valve located between IHX and the 
reactor vessel limits the sodium flow rate during decay 
heat removal operation, thus reducing thermal shocks. 

All pipes except those in the annular region surr- 
ounding the reactor vessel are provided with an electric 
trace-heating system. The pipes in the reactor vessel, 
the pumps and iHX^s are preheated by hot nitrogen. 

The pipes of the three primary loops are laid out 
symmetrically in the reactor cell, from which they run to 
the three parallel primary cells. This arrangement cases 
shielding of the piping penetrations and reduces the 
activation of the structural material and of the secondary 
sodium. 
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2,2 Fault Tree Construction 

First of all referring to 3^R~300 figure we divide 
loop into two parts, 

(a) Primary Loop - Consisting of pipe segments 
1,2, 3,4, primary pump. Butterfly valve and heat exchanger, 

(b) Secondary Loop - Consisting of pipe segments 
5,6,7, shut off valve, superheater, evaporator and secon- 
dary pump. 

Then following observations were made* 

(i) Pump has two kinds of failures - (a) single 
ended or minor failure? (b) double ended or major failure, 

(ii) Different pipe segments can in general (and 

in fact they are so) bo of different diameters hence their 
failure data is expected to be different from one another. 
As a result, these pipe segments break or leak must be 
considered as separate basic events, 

(ill) Heat exchangers. Butterfly valves, primary 
pump, secondary pump, shut off valve, supeiiieater, 
evaporator have redundant duplications intended to improve 
reliability, 

(iv) Only when primary or intermediate heat 
transport system and emergency core cooling system fail, 
there is uncontrolled rise in core temperature. 
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(v) Since major components have redundancies, 
there must be maintenance vtien one or more of redundant 
and main component fail. Hence repair rate must be 
taken , 

(vi) Primary pump is a sodium pump and is more 
sophisticated. Its data must be chosoa carefully., 

(vii) There are three identical primary loops hence 
over all system is 1 out of 3, 

Based on above information we make further sub- 
division of the system as follows; Primary loop is 
divided into two parts, (a) Valve line consisting of pipe 
segments 3,4 and Butterfly valve, and (b) Pump line 
consisting of pipe segments 1,2, primary pump and inter- 
mediate heat exchangers. Any of these parts’ failure will 
lead to primary loop failure. Valve line failure occurs 
when either there is no supply from Butterfly valve or 
pipe length 4 has failed. Butterfly valve supply is 
cut off when either pipe length 3 is cut off or Butterfly 
valve IS stuck. Pump line failure occurs when either there 
IS no input to heat exchanger or heat exchanger fails. 
There is no input to heat exchanger if either pipe segment 
2 fails u or no supply from primary pump. There is no 
supply from primary pump if pipe segment 1 fails or pump 
fails, 

Inteimediate loop fails if either valve line 
(consisting of pipe length 6 and 7, shut off valve and 
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sec, pump) fails or super heater, evaporator line fails 
(consisting of superheater, evaporator or pipe length 5), 
Valve line fails if either pxpe length 6 fails or sub- 
system from valve to heat exchanger fails. The later fails 
if pipe length 7 fails or pump and valve systan fails. 

Pump and valve system fails if either of these fails. 
Superheater, evaporator line fails if pipe length 5 fails 
or any of superheater and evaporator fails. 

Emergency core cooling system fails if either air 
cooler fails or there is no supply to air cooler. There is 
no supply to air cooler if either emergency pump fails or 
pipe length fails. 


2,3 Data to be Used 

Based on Appendix I, we decided to use the following 

data; 


Component 

X hr”^ 

hr"^ 

Emer. cooling pipe (single ended 
failure) 

2.0x1 

1 .2x10”^ 

Emer, cooling pipe (double ended 
failure) 

2.0x10^'*° 

1 .0x10”® 

Other pipes (single ended) 

1 .0x10'"''° 

5.0x10”^ 

Other pipes (double ended) 

1,0x1 0 '-^^ 

4,0x10”^^ 

Valve (Butterfly) 

1 ,6x10“^ 

1 .0x10”"^ 

Primary pump* 

O.ITxlO"*'^ 

1 .0x10*"^ 

Heat exchanger 

4,184x10”° 

2,0x10”^ 
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Compon ent 

K hr 

hr~'' 

Secondary pump 

2.5x10''^ 

2.5x1 0"^ 

Shut off valve 

3,0x10“^ 

3,0x1 

Super-heater 

8,0x1 

5.0x1 0“^ 

Evaporator 

8,0x1 CT^ 

5.0x10'’^ 

ECC pump 

2.5x1 0~^ 

2.5x10“^ 

Air cooler 

3,0x10“^ 

5.0x10“^ 


2,4 Reliability Program Objectives 

We aim at ensuring that the liklihood of exceeding 
the nuclear radiation dose guidelines of 10CFR 100[ ] 

at the plant site boundary should not be greater than 10 ° 
per reactor operating year ^hat in-place core-coolable 
geometry will be lost [ & ] , Loss of in-place coolable 
core geometry is a failure criterion which is used to 
characterize very lov; probability events. If core coolable 
geometry is lost, there is no assurance that significant 
core damage could not occur, even though there is still a 
lav probability that site boundary dose guidelines could be 
exceeded S3nce the plant includes a number of containing 
barriers, one of which is the containment structure. The 
probability of occurrence of potential initiators of loss 
of coolable geometry can be controlled by design. It then 
becomes the task of the Reliability Program to assure the 
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hi^ reliability of systems necessary to prevent the 
onset of such initiating events. Basic tasksof the 
Reliability Program are (1) to identify those extremely 
unlikely events having the potential to produce loss of 
coolable core geometry (2) to ensure through reliability 
engineering design that all such events are of sufficiently 
low probability to meet the goal set (3) to confirm this 
high reliability through amlytic assessment and testing. 

An initial reliability allocation against the 
overall goal was made for potential initnators which might 
threaten coolable geometry [ 8 ] , The division of the 
overall numeric goal among safety systems was based on 
the functional role of each system and its predicted 
relative failure potential. On the basis of this, the 
preliminary unreliability allocation established for shut- 
down heat removal system [BJ was < 8x10 failures per 


year 






















CHAPTER 3 


PROGRAM PCOMCP FOR FINDING MINB4AL CUTSETS 


3.1 Definition 

Cutsets— A cutset is a colloction of basic eventsi if all 
these basic events occur, the top event is guaranteed to 
occur. 

Pathset - It IS a collection of basic events and 
if none of the events in the set occur, the top event is 
guaranteed to not occur. 

Minimal cutset - A minimal cutset is such that if 
any basic event is removed from the set, the remaining 
events collectively are no longer a cutset. This is 
obtained from cutset by removing cuts which are redundant 
according to some rules laid dovn below. 

Minimal Pathset - is a path set such that if any 
basic event is removed from the set, the remaining events 
collectively are no longer a pa-th set, 

3.2 Redundancies to be removed for Generation of ivlinimal 

tiuts'et ' ' * ' ' ^ ^ . ,, , „ 

(i) Redundant Factor 

If in any cut we have A. A, A or A like term 


replace it by A 
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(ii) Subset Redundancy 

If two cuts are such that one is subset of other 
remove the superset, i,e, A + AB = A, 

(ill) Term Redundancy 

If two terms are repeated in the cutset then 
remove the repeated term, AB + AB = AB, 

3,3 Algorithm 

The program is based on the unique factorization 
properties of natural numbers stated in the form of 
following theorem. 

Unique Factorization Theorom f 1 

Every natural number greater than 1 can be expressed 
as a product of prime factors in one and only^way, apart 
from the order in which the factors are written. 

Our strategy is to assign one prime number to 
each basic event, doing arithmetic operatioriBBfor logical 
GATES, removing redundancies in cutsets and finally decode 
them back, 

A particular combination of basic events can be 
expressed uniquely as a single number whidi is equal to 
the product of prime numbers corresponding to the basic 
events. 



15 


3,4 Program Details 

We follow the follov\/mg steps? 

(i) Assignment of one prime number to each basic 

eve nt , 

(ii) Go bottom up, 

(ill) Go up level by level (to be described), A level 
IS the number of gates from top after which the event is 
encounte red, 

(iv) For each of OR gate encountered we preserve 
all the inputs to gate in an array, 

(v) For each AND gate we multiply the input codes, 

(vi) Of the final array, if two number is multiple of 
or equal 

other larger^is removed. This is to remove of redundancy 
No, (ii) and ( in) , 

(vii) Of the final set, eacii numbers is prime 
factorized. If any factor isrepeated in the factor set 
they are removed. The factor set when decoded gives minimal 
cutset. This step removes the redundancy (i), 

(viii) If number of prime factors increase the 
maximum order of cutsets wanted, that cut is als'o removed, 

(ix) AND gate is designated by an even number to 
the gate and OR gate by an odd number. 

Method of inputting the tree structure is by 
assigning a series of integers w^ich represent the gates 
linking the top event to a primary event. The length of 
each senes is equal to maximum level (i,e, maximum number 
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of gates from Top event to basic event) plus one. The 
gate number 0 indicates that there is no gates at that 
position. The order of row corresponds to the ctnfigura- 
tion of a fault tree. 

Level by level is done as follows' 

First line of input data is read and level is 
decided. The number is stored in the matrix element 
N(L,1) where L is level number. Second line of input 
data is thai read. If level number is same as previous 
number, the last gateo operation is done. If last gate is 
OR gate (i,e, gate number is odd), it is stored in N(L,2) 

1 ,-e, in the N(L,x) which is vacant (i.e. zero) and x is 

minimum. On the other hand if gate is AND gate (i,e, gate 
number is even) all N(L,x) elements are multiplied by 
the current line basic event number. Once the gate 
operation is done level is reduced by one and this number 
IS stored. Now the third line is read. If the level 
number equals tho level number stored, the gate operation 
is done according to rules aforesaid and gate type 
indicated by the gate number being odd or even at the level’s 
place from the Top gate number row wise. Level number is 

reduced by one and this level is stored. In cases when the 

current level is not equal to the stored level number, the 
current basic event coded number is stored in N(L’,1) and 
next line read, if level is again not equal to previous 
level number, this basic even number is stored in N(L”,1) 
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and this is carried till the two successive levels are 
equal. Gate operation is done and level number is again 
reduced. If this level pointer is equal to previous level 
number, proper operation done and this process is repeated 
tiJ 1 all the numbers are in N(1,x) where x varies from 
1 to No, of OR gates in the fault tree. 

To remove the redundancies (ii) and (iii), the 
final set of numbers N(1,x) V x <No, of OR gates, is 
checked for multiplicity. If any number in this set 
IS multiple of any other, the larger (or equal) is set 
equal to zero. Once "this procedure is carried out 
supersets are removed. Next for eadi non-zero elements 
in N( 1 , x) set, prime factors are found. If any prime 
factor IS repeated in any element in N(1,x), this repeated 
factor IS removed, i.e, one of repeated factors are 
kept. This process ensures the removal of type (i) redun- 
dancy, Corresponding to each number in setN(1,x), the set 
of non-repeated prime factors when decoded give minimal 
cutset elements. The number of minimal cutsets equal the 
number of nonzero element in set N( 1 , x) , 

3.5 Salient Features of the Program 

(i) Coding by prime numbers ease out the minimi- 
zation of cutset procedure, 

(ii) Logical operations carried out by anthmatic 


and matrix operation 
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(iii^ Actual basic event No, when coded by prime 
numbers does not increase the CBE value very much as 
compared to actual basic event Ho, ( CBE=Coded Basic Event). 

(iv) The maximum order of cutsets could be specified, 

(v) Level by level operation is actual practice 
reduces the total number of operation required, 

(vi) Because of (i)| (m) and (v) the program 
IS efficient and takes loss CPU time, 

(vii) Memory storage required is loss than many 
other minimal cutset enumeration programs, 

(viii) The maximum number of OR gates and input right 
now could be 100, This can bo increased by simply increasing 
the dimension of some arrays in the program, 

3 ,6 Limitations of the Program 

(i) All gates must have only two inputs. Only 
the top gate can have 3 or more inputs, 

(ii) As it IS the program takes care of only AND 
and OR gates. Some more gates can be taken into account 
as mentioned belows 

(a) NOT gate - the basic event itself is changed. 

It IS considered independent event vdth R = 1 ~ not inverted 
event reliability, 

(b) INHIBIT Gate - It is considered an AND gate. 

With conditional event as a basic event. 
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( c) Priority AND It is bosically AND gate with 
some sequoncingj so it is treated as AND gate, 

3 ,7 Operations to be Carried Out 

The operation can be summarized as below; 

(a) A line is read and level is decided, 

(b) If level IS equal to previous step level no. 
Logical operation is to be done otherwise CBE is stored in 
N(L,1), 

(c) Logical operation is carried out as follows; If 
it IS an OR gate the CBE is to be stored in vacant place of 
N(L,x). On the other hand, if it is an AND gate all the 
elements of N(L, x) are multiplied by CBE, 

(d) If any logical operation is done LP is decreased 
by 1 , If NG (LP) is even,(i,e, AND gate) all combination of 
multiplication of N(LP, x) and N(LP +1, x) is carried out 
and stored in N(LP, x). In cases when NG (LP) is odd, 

All elements of N(LP + 1,x) stored in vacant (i,e, zero) 
places of N(LP, x) , All elements of N(LP + 1, x) reduced 
to zero, 

(e) LP decreased further and ( d) line operation 
carried out till LP is zero, 

(f) When all of input data is over, factorization, 
minimization and decoding is done as outlined above. If 
data is not over go to (a). 



CHAPTER 4 


SYSTEM QUANTIFICATION ~ A THEORETICAL OVERVIEW 

4.1 Some Definitions F »3] and Symbol Meanings 

(a) Reliability at time t = R(t). 

The probability that the component experiences no 
failure during the time interval (O,"!] given that the compo- 
nent was repaired at time zero. 

(b) Unreliability at time t = F(t) = 1-RCt) (4 

(c) Failure density at time t = f( t) 

The first order derivative of F(t). 

( d) Failure rate = r(t) , 

The probability that the componoit experiences a failure per 
unit time at time t given that the ccxnponent was repaired at 
time zero and has survived to time t. 

(e) Mean time to failure = MTTF 

The expected value of time to failure, i.e, mean 

of the span of time from repair to first failure 

00 

MTTF = / t f(t) dt (4 

o 

(f) Repair probability at time t = G(t), 

The probability that the repair is completed before 
time t given that the ccxnponent failed at time zero, 

(g) Repair density of G( t) = g(t) 

The first order derivative of G(t). 
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(h) Mean time to repair = MTTR, 

The expected value of time to repair, i,g, mean of 
length of time from the failure to the succeeding first 
repair 

00 

MTTR = S t g(t) dt (4.3) 

o 

(i) Repair rate = m( t) , 

The probability that the component is repaired per 
unit time at time t given that the component failed at 
time zero and has been failed to time t, 

(3) Availability at time t = A( t) . 

The probability that the component is normal at time 
t given that it was as good as new at time zero. 

A(t) R(t) . 

( k) Unavailability at time t = Q( t) 

The probability that a component is in the failed 
state at time t, given that it jumped into the normal state 
at time zero 

Q(t) < F(t) 

(l) Conditional failure intensity = ( t) 

The probability that the component fails per unit 
time at time t, given tnat it is in the normal state at 
time zero and is normal at time t. 

(m) Unconditional failure intensity at time t = w(t) 

The probability that a component fails per unit time at 
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t, given that it jumped into the normal state at time zero, 
(n) Expected number of failures ( HIF ) = W(t,t+dt) 

Expected number of failures during [t, t+dt) , given 
that the component jumped into the normal state at time zero, 
(oj W ( t^ , = ENF over a period. 

Expected number of failures during ( t^ , i given 
that the component jumped into the normal state at time zero, 

WCt^yt^) = /'^w(t) dt ('^.4) 

(p) Conditional repair intensity = P'(t) 

The probability that a component is repaired per unit 
time at time t, given that it jumped into the normal state 
time zero and is failed at time t, 

( q) Unconditional repair intensity at time t = v(t) 

The probability that the component is repaired per 

unit time at time t, given that it jumped into the normal 
state at time zero, 

( r) Expected number of repairs is an interval 
= V(t^, t2) 

Expected number of repairs during C > "^2 ^ ^ given 
that the component jumped into the normal state at time 
zero. 

"^2 

= / v( t) dt 

t. 


VCt^ft^) 
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4.2 Fundamental Relations f ’^ 1 

A(t) = R(t) for non-repairable components 
X(t) = r(t) for non-repairable components 


R(t) = / F(uj du 

t 


r(t) = 


-Fit. 


(4.6) 


(4.7) 


R(t) 


m( t) = 


exp [- £ r ( u) du ] 

o 


f(t) = r(t) exp [- / r(u) du ] = r(t)R(t) 


q( t) 
i ~g( t ) 


G(t) = 1 - exp [“ 1 m(u)du] 


(4.8) 

(4.9) 

(4.10) 

(4.11) 


g(t) = m( t)[l-G( t) ] = m(t) exp[- ■/' m(u)du] (4,12) 


If constant failure rate =^= r(t) and non-repair able 


X ( t) = 


~R( tj 


(4.13) 


R( t) = e' 


(4.14) 


F(t) 

f(t) 


1 - e 
. -H 


(4.15) 

(4.16) 


MTTF 


(4.17) 


If constant repair rate 
G(t) = 1 - 

g(t) = \i 


p- = m(t) 


(4.18) 

(4.19) 


MTTR = 


(4.20) 
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4,3 Relations among the Whole Process 

(a) Unconditional intensixies w( t) and v( t) . The 
component which fail during (t, t+dt} are of two types; 

Type 1 ;A Component which was repaired during 
[ u, u+du ) , has been normal at time t, and fails during 

[ t, t + dt ) , given that the component jumped into the normal 
state at time zero. 

Probability for such component is v(t) du,f(t-u)dt 
since v(u)du)= probability that the component is repaired 
during [u, u+du), given that it is as good as new at time 
zero and f(t~u)dt = the probability that the component has 
been noimal to tine t and failed during ( t, t+dt] , given 
that it was as good as new at time zero and v/as repaired at 
time u. 

Type 2 ; A component v\hich has been normal to 
time t and fails during [ t, t+dt) , given that it jumped 
into the normal state at time zero. 

Probability of second type of components is f(t)dt 

w( t) dt = Probability that the component fails during [t,t-^dt ) , 

given that it jumped into the normal state at time zero 

t 

w(t)dt = f(t) dt + dt f(t-u) v( u) du 

t ° 

or w(t) = f(t) + f( t-u) v(u)du (4,21) 

(b) Second Relationship; The component which ore 
repaired during [t, t+dt) is one which was failed during 
[u,u+du) I has been failedto time t and repaired during 
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[t, t+dt ] , given that the component jumped into the normal 
state at time zero. 

The probability of such components is w( t) du,g(t ~u)dt 

t 

v( t) dt = dt / g( t-u) w(u) du 
0 

* t 

. . v( t) = / g(t-a) w(u) du (4,22) 

0 


(c) Unavailability? 

q(t) = W( 0 ,t) - V( 0 ,t) (4.23) 

= Number of failures - Number of repairs 

at time t 

( d) Failure intensity ^(t) , 




(e) Repair intensity p.(t) 

nCt) = ^ (4.25) 


4 ,4 Laplace Transform Analysis for the Whole Process 
Laplace Transform of (4,16) 

L[f(t)] = 


S+X 


L[g(t)] = -g" (4.19) 

Laplace transform of (4,21) and (4,22) 

L[w(t)] = L[f(t)] + L[f(t) ], L [v(t)] 

L[v(t)l = L[g(t)] , L[w(t)] 



26 


Using the L [ f( t) ] and L[g(t)] 

L[w(t3 = ^ [v(t)] 

L[v(t)] = LCw(t)] 

Solving the two equations for L w( t) and L v(t) 

L[w{t)] = ^ (^) + ^ 

L[v(t)] = ^ (i) - (s 5 ^) 

Taking Laplace inverse 


w( t) = 

Xp, 

X+p 

-f 

x^ 

x+p 



(4,26) 

y( t) * 

Xp 


xp 

x+p 

e~(X+P)t 


(4,27) 

W(0,t) = 

Xp 

X+p 

t 

+ ■ 

(X+p)^ 

(l-e‘ 

.(X+p)t_^ 

(4.28) 

v(o,t) = 

Xp 

x+p 



Xp 

(x -tp 

(1-6 

3,“(X+p)t J 

(4.29) 

Q( t) = 

w(0,' 

t) 

- v( 0 ,t) = 

_ 2 ^ 
?v + u 


) 








(4.30) 


4.5 System Analysis 

A, Component analysis is done on the basis of 


following flow chart; 



i(4.26,4.27) 



Unconditioned intensities 
w( t),v( t) 


(4.4,4,5) 



Modified version of above flow chart will be discussed in 
the next chapter v\tiich we used in our program, 

KITT - The code is an application of kinetic Tree theory 
and handles independent basic events wnich are non~repairable 
or repairable, provided they have constant failure fates 
and constant repair rates 
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B, Flowchart for KITT Compute 



Exact time-dependGnt reliability parameters are 
determined for each basic events and cutsetsv-~t>ut for the 
system as a whole the parameters are obtained by upper or 


lower bound approximations or by bracketing. In bracketing 
procedure, the various upper and lower bounds can bo obtained 
as close to each otiier as desired, and thus the exact value 
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for system parameters are obtained if the user so choses, 

w(t), v(t) are obtained by equations (^,21 ) and (4,22), 
For this purpose, numerical integration is used, W(0,t), 
V(0,t) are obtained by another integration for given t, Q is 
found by equation (4,23), 

Minimal Cutset Parameters? 


(a) Unavailability— A cutset is occurring if all 

the basic events in the cutset are occurring 

n 

Q * (t) = 71 Q.(t) (4,30a) 

3=1 ^ 

(b) Conditions failure intensity and Unconditional 
failure intensity; 


X * (t) 


w * ( t) 
1- Q*(t) 


(4.31) 


where 


w * ( t) 
(c) Similarly 

V * (t) 

* (t) 


n 

E w(t) 

1=1 ^ 


n 

It Q. ( t) 
1=1 ^ 
15^1 


^ v(t) 
1=1 ^ 


n 

7C 

1=1 

1?^1 


[1-Ql(t)] 



( d) Also 

W * (0,t) 
V * (0,t) 


t 

/ w * ( u) du 
o t 

/ V * (u) du 


(4.32) 


(4.33) 

(4.34) 


(4.35) 

(4.36) 
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System Parameters * 

(a) Unavailability 

Let = event that all the basic events of ith 
minimal cutset exist at time t. 

The ith minimal cutset failure exist 
at time t. 

By inclusion-exclusion principle 

^c 

"t) = Pr ( U cL ) ; N = No, of minimal 

1=1 cutsets 


n N^ 1-1 

C C 

= S Pr( d J - S 2 

1=1 ^ ±=2 3 -- 


\m 


+ ••..+ (-1) E 

1 

N -1 

+ •••• + (~1) ^ 


Pr(d.(\d ( + 

1 ^ ^ 

PT(d^,nd,2n-..H„) 
i ''c 

Pr(d^f\ d^n ) 


1 < < i2.^,i < N 

* • • Tr » X — 


The mth term is the contribution to Q_(t) from m minimal 

cut set failures existing simultaneously at time t 
N_ 


Q (t) = S Q * (t) - 2 

s 1=1 ^ 


Nc 1-1 

2 ,^Q(t)+.... 

1=2 j =1 



, 2 11 ."I. 

N -1 


+ (-1) 

/ / Q( t) 


m 


(4.36a) 


where 


l '^ • • ol ; 


IS the product of Q(t)‘s for the basic 


m 


events in cutset i^ or io.«...or 


'm* 
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Bracketings If in above expression, only first teim is 
taken, it is upper bound for Q (t). If first two terms 
taken, it is lower bound for Q^Ct), If 1st three terms 
taken, it is better upper bound for Qg( t) and so on. 

Upper bound alternative expression is due to Esary 

N. 

and Proschan = 1 - % [1 - Q *(t)] , 

1=1 

This upper bound estimate is sometimes conservative 
estimate, but it is exact when the cutsets are disjoint sets 
of basic events, 

(b) Unconditional failure intensity 

t) = Expected number of times the top event 
occurs at time t, per unit time 

w(t;1,.,,m) = The uncondational failure intensity 

for a mode of failure v;hich has as its basic 

failures the basic failures which are common 

members to all the mode failures 1,..,fm 

/ . \ N 1-1 

(t) = w *(t) - 2^ 2 w(tji,j) ^ Q(t) 

s 1=1 ^ 1=2 j =1 

1-1 J-1 ^ Ti; , V 

+ 2 2 w(tj i»J,k) , T V t) + ,,,, 

1=3 J=2 k=1 

N -1 

+ (-1) ^ w(t;1,...,N^) Q(t) (4.38) 

A = The event that one or more of cutset 
failure occur at time t, 

B = A rv X 

WgC tj1, dt = Pr(e^ f\ f\ej^n B) 

H N 1-1 

= Pr(e.f\,..Oo^^r\d^) - 2 ^" 2 Pr( e^O. . .ae^^n 

1=1 ^ 1=2 3=1 

_nd Hd (-1) ° Pr(e Qd ) 

J. J I Q 
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w 


( 2 ) 


Nc x-1 

2 w(tji)dt- 2 2 w(t;i,j) + . 

1=1 ^ 1=2 j =1 ^ 


N -1 

+ (-1 ) 


WgCt) 


Bounds; 


( l) W ( t) 

s min 

= ^ ^ ( t) - w^^^(t) 

s ^min s ^max 

w^( t) 

'^max 

= ^ t) ~ -h t) 

( ii) if N la 

Vi# 

even* 

Q,(t) = 

Q (t) 

s^ ^min 

w^(t) = 

'^s^ ^^min 

A3(t) = 

^s^ ^^min 

If N IS odd; 
c 


Q^Ct) = 

Q (t) 
s ^max 

W3(t) = 

w ( t) 
s'' ^max 

\s(t) = 

Xi _ V /y,n ^l. - . 

s max 

(ill) 

w 

= _^max 

1 

^Sfmax 


X ( t) 
s ‘^min 


w 


s»min 




(c) Integrated number of failures 

t 


Wg(0,t) 


/ w ( u) du 
o ® 


(4.40) 


(4,40a) 

(4.40b) 


(4.41) 
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4 .7 Alternative Formulation of System Analysis ~ Short-Cut 
Technique 

J,B, Fussell formulated this technique .which is 
basically back-of the~envelope guesstimate. It requires 
as input failure and repair rates for basic events and 
minimal cutsets. It assumes exponential distribution of 
and and independent of component failures. 


(a) Component Level Analysis ; 

-Xit 

= 1 - e ^ = X^t 

if X^t < 0,1 
If components repairable 


Q 

^1 




)t 

(1 - e ^ ) 


as t becomes large and if 

Q rJ 


^1 X^ 

— i << 0,1 

^^1 


A. 

1 


1 




i if t > 


t4,42) 


(4.43) 


( b) Cutset level Analysis ; 



K t) = 

n 

1=1 ^ 



(4,44) 

Using equations 

(4.32) 

and (4,24) 

V 





n 


n 


'^1 * 

(t) = 

S i[1 - 

j=1 

Qj(t)] 'X3_(t) 

% 

1=1 

Qj(t) 








Substituting equation (4, 44) here, 
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or 


n 

* (t) ~ * (t) 


noting 1 -* Q^(t) « 1 


K 


* (t) = 


1 - *C^) 


■^(t) 


(4.45) 


(4.46) 


( c) System Level Analysis using bounding Procedures ; 


Qjt)co 

s' * ( t) 

(4.47) 


1=1 





X CV» 

s — 

s'" X ^ 

i=1 1 

(4,48) 

W (■fc).o 

s **• 

Nr. 

s ° w * ( t) 

1=1 

(4.49) 

s 

W3( t) 

T^^t) • 

(4.50) 



CHAPTER 5 


PROGRAJVt EOR SYSTEM QUAMTIFICATJON t MODKITT 


5.1 Formulation for m Out of n Identical Systems 
(a) Unavailability’s 

m out of n s’ystsm: Top event unavailable if any 
of m or more systems out of n systems are unavailable. 

By binomial dis tnbu’tion. 


Q,(t) =2 ( k) ^ 

® k=m 


(5.1) 


(b) Unconditional and conditional intensities for m 

n n 

out of n s'ystems” There are mth order cutsets; (jn 4 .-|^» 

(m+1)th order cutsets and so on. If we noglect higher order 

n 

cutsets we are left with mth order cutsets. For finding 

miniiml cutsets, we must do so. 

Using equation (4,32) for any cutset of mth order 

w *( t) = m w 

ni“" i 

Since there will be m identical -teims with valve w Q 
Finally using Fussell’ s approximation [‘3] equation (4,49) 


W^(t) 




I ^ni-1 

*""-r m Q w 




ml nHoil 


Ws(t) 




nl rvin-l 

. Q 

m-1 I n-ml ^ 

Wg(t) 


(5.2) 


(5.3) 
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(c) Unconditional and conditional repair intensities 
using Gqn.(4,33) 


V * (t) = rn v( t).( for any of 


cutsets of mth order. 

Using Fussell' s approximation 


v^Ct) 


{„) m V (1-0)‘ 


e-.Ct) 


m-U b~m'. 

V { t) 


(5,4) 


(5,5) 


5,2 Flow Chart for Svstcan Analysis 

Before looking at final flow charge, let’s have a look 
on the ccmponent or basic event analysis. 


Data 



(4,4, 4, 5) (numerical 

integration) 


Expected numbers 
W(0,t), V(0,t) 


(4,30) (Numerical 
iteration) 

\/ 


Conditional intensity 
p- ( t) 
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5 ,3 A Note on m Out of n System 

If we do not define m out of n systan as done in 
section 5,1 rather as ”Top ©vent available if any 
of m or more systans out of n systems are available" - 
we can still use the formulation of (5,1 ) if we do the 
following. Replace m by nnn+l , 









^ -^w I ^ CWfjsntvsl v=\, ,t» 



•’ >-W)'l) , \lo,\) , 61 i,Ot) J > ^4), 

«^W/>A^ ivtmrx -rta. <Nbl>\t C^Urt 


V 

CAAjfc-^ ® 

V4-H (^J 


t4- 3o«t) 


Ciuhs^t 
Q.-* Ct) 




1(4 ‘SO 
CjudcSA^ 

(3^ 


1 

6.<, T^^ecw 


1(4- SO) 


S^Wx ‘\<ar>l 


it -Ha. 'Vj*» 

TY\ ov;;*: O-lj TTN S>^V'<AV\^ 

<j\)^t/Tvsrosji_ M"Crb 


TrW-X y HVN-VM ^ VVNftjCTVI 

XcB ij •> m ‘ “ 


^ tfe 
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© 



Not all cutsets are taken. There are 3 optionss 
Option 0 => Take all cutsets 


Option 1 =^> Maximum order of cutsets to be considered 
are specified, Hii^er order cutset Qts^ 
are set to zero. 


Option 2 Or more => /^ accuracy desired is specified, 

0* - Q * 

If — > % assigned 

max 

Q * IS set to zero. 


Bracketting has been used to find Q_ _ » 

6 f ms X 

I^" "the first three terms taken 

S f III J.I i 

in equation (4,36a), it gives ,^3^* 
first two terms taken, eqn,(4,36a) gives 
Q 

^s , min • 


Q mean 
s 


V-Q * 0 

s,max s,min 


(5,6) 



This geometric mean has been suggested in 
WASi--1400 report [4,5] , 

For all these quantities 

Qo same equation (5,1) has been used, 

S y UlSali 

Hence and are not very accurate 

o j lUa X SyluXXi 

but Q IS exact, 

^s,mean • 
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5 .4 Input to Proarnm 

In file nomed FOR 2^. DAT gives 

In the first line — Time ?ind decimal accuracy of 
result wanted. 

In 2nd line onwards component no, ^(Q), M,N, 

MN for all basic events, M,M indicates M out of N system, 

M = 0 => PI = X P2 = p 

MN ^ 0 => PI = Q(T) P2 = w(T) 

Next line must begin with 0 to indicate end of basic 
events , 

Next line we have to give OPTION = ? 

If Option = 0 => take all cutsets. 

If option = 1 => Next line must contain 

MAX ORDER = ? 

If option ^ 2 => Next line must contain % = 7 

Next lines contain cutsets one in each line till all cutsets 
are over. End of cutset information again indicated by a 
line starting with 0 or a black line. 

Last lino has M and N for the TOP GATE, 

5.5 Output 

First few lines have w(t), v(t), W(0,t), V(0, t), 
Lambda, Meu. and Q( t) for basic events. 

Next wo have 

w * (t), X * (t) and Q * ( t) for the relevant 


cutsets only 
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Lastly, w(s), Q(s)min, Q(s)max, Q(s)mean for the total 

sys tem. 

Note; W(0,t), V(0,t) is set to zero for some basic events, 
this has been done to reduce computational effort. These 
quantities are required for calculation of Q(t), so i f 
Q( t ) IS obtained by any alternate formulation or as such 
provided in input these quantities have no relevance and 
hence have been set to zero, 

^ Merits of the Program 

(i) In all formulation the most important quantity 
Q(t) (unavailability) calculation has been done exactly except 
for bracketing. First three terms of equation (4,36a) has 
been token for upper bound wh3ch make it mors close to actual 
value. Even for M out of N system unavailability calculation 
IS exact, 

(ii) Fussoll's approximation has been used for quanti- 
ties which are not of primary interest like w(t) of m out of 
n system and w (t) for the entire system. Rest of the 

o 

other quantities* formulation is exact likeW(0,t), V(0,t), 
w(t), v(t) for componaits. 

(iia) Answer accurate to any decimal place can be 
obtained. The iteration procedure and integration procedure 
has been written in a such a manner that it automatically 
reduces the step size of iteration for obtaining the desired 
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decimal place accuracy, 

(iv) Only those cutsets hove been considered whose 
reliability is large enough according to certain criterion 
enumerated in the form of OPTION, This reduces the compu- 
tational effort drastically. This reduces the computational 
effort drastically. In the current input, onlv 12 out of 
total 84 cutsets are important for 99% accuracy, 

5 ,7 Demerits and Suggestions for Further Improvement 

(i) The bracketting procedure has not been kept flexible. 
The number of terms for upper bound approximation and lower 
bound approximation has not been kept variables but aro 3 

and 2 in this program. So desired accuracy of not 

be specified. This is not a particularly severe demerit 
since the successive terms of equation (4,36a) decrease by 
Q( t) which IS a small number and hence two bounds found above 
IS almost close to the actual value . In any case is 

found by a method suggested by WASH - 1400 report,H9nce, it is 
very dependable, 

(ii) Fussell's approximation has been used for some 
quantities which is valid only for large time. This again 
IS not a severe constraint since the approximation has been 
used only for quantities of secondary importance. MTTF 
and MTTR are in a few hours generally less than 100 whereas 
the time to be considered for calculation is in years. Hence 
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condition of the kind equation (4,43) and large time is 
usually satisfied, 

(ill) v(t),p. (t) for cutsets and overall system has 
not been found. This has not been done in KITT either. 

The reason for this is that these quantities are not of 
primary interest. Secondly, not all components are repair- 
able hence calculation in those cases is not necessary, 

5 ,8 Discussion of Result 

For LMFBR primary loop, we obtained the following 
result for unavailability of the svstem fault tree we 


started 

out with , 




- 

Mjl 

A ( s) 

w( s) 

Gl . s )min 

Q(s) 

max 

Q(s) 

mean 

1 year 

0.b916x10~^^ 

0,5916x1 0”"'^ 

0.25666x10"® 

0.257X 

10"® 

0.2569X 

10"® 


. 0.9374x10"''° 

0.9374x1 0"''° 

0.7857x10"'^ 

0.7938X 

lo-'^ 

0.7898X 

10""^ 


After some improvement in the system by increasing redundancies 
in components (in parallel) as in the pnnt-out, we got the 
followings 
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Time 

WBt 

w( s) 

Q( s )min 

Q(s) 

max 

Q(s)“ 

mean 

1 year 


0,8637x10"^^ 

—1 0 

0.38851x10 

0 ^38868x 

10 

0.38859X 

10-1'^ 

40 years 

! 

j 

i 

0,3425x1 O"''^ 

0.3425x10"^° 

1 . ^ 

0,24798x10'’^° 

0.25016X 

lo-'^ 

0,24907> 

10"’’^ 


Wg observe the following facts from our analysis s 

(i) With time, system unavailability increases , Even 

—7 

at the end of 40 years unavailaoility is 0,25x10 to 0,8x10 
which is quite satisfactory, 

(ii) As we increase redundancies, unavailability gets 
reduced. Our analysis can, thus, help in system design. We 
can decide the redundancies (standby and parallel 

components) based on some top event failure probability 
according to certain unavailability allocation techniques. 
Adequacy of the design can also be ascertained, 

(ill) Even for 99?^ accuracy, not all minimum cutsets 
are important. So, it is better to concentrate on the 
important minimum cutsets rather than all. Also with time 
the important minimum cutsets keep on changing. Some old 
important cutsets no longer remain that important and seme 
new adds to the list of important ones. 


(iv) As the % accuracy required is decreased, not 
many new cutsets are added. So we can assign less accuracy 
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desired without much increase in error. Like, in place of 
99% t we can assign 85/^, 

(v) Some components are more critical in determining 
the system unavailability. Change of redundancy of sane 
component has much pronounced effect on the overall unavail- 
ability than others. 

-7 

Our result is closely comparable to the value 4x\0 

per roactor year obtained by F.J. Baloh, N,'/V. Brown, J, Graham, 

A.M, Smith, P.P, Zemanick [S] for Clinch River Breeder 

reactor plant. The value is also consistent with the primary 

-7 

allocation (goal) set which was loss than equal to 8x10 , 

This number was allocated on the basis of 10 CFR 100 criterion 
near plant site boundary [ S ] , 



CHAPTER 6 


ILLJSTRATIVE PROGRAM FOR MINIMAL PATH AND CUTSET 
FOR 'AN 'ELECYrICAL sV^'BA 


6,1 Introduction 

Electrical systems differ from the basic nuclear 

system in that its fault tree has a network graph like 

structure. There are many feedback paths, parallel paths 

and interconnections. This is not like the fault tree 
containing OR and AT^D gate and hence a separate formulation 

IS necessary for its analysis. 

In the technique usea in the program, we do not 
require all the minimal paths to be deduced and checked, 

A few paths (called basic minimal paths) are deduced from 
the minimal path tree. The combination of failures that 
breaks the set of basic minimal paths is sufficient to 
deduce all the minimal cutsets of the n/w. The set of 
basic minimal paths is a subset of all the minimal paths 
of the n/w and the remaining minimal paths are not necessary 
for evaluating the minimal cutsets [ao], 

6,2 Definition 

Minimal Path - It is a path from source to sink 


whereby no nodes are traversed more than once [ 13 ]. 
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Basic Minimal Path — It is a minimal path in v\tiich 
no element of the path are linked to another by any branch 
in the n/w except with ihose elements before or after it 
in the path, 

6.3 Program Details 

Program is closely based on [ algorithm details of 
which can be looked into , What we have done is an 

improvement in it to make it more efficient and doing min»r 
corrections. The output formal has bean changed to make it 
more useful. 



CHAPTER 7 


FAULT TREE MODIFICATION PROGRAM (PR0GRAM~4) 

7 .1 Ob.iect 

Apart from AND and OR gates, there are many other 
gates which occur in the fault trees like XOR, NOT, NAND 
NOR, Priority AND, Inhibit etc. Priority AND and Inhibit 
gates are treated as AND gates. There should be some method 
to take care of XOR, NOT, NAND, NOR, M out of N gate, 
out of *N* gate we took care at quantification level. The 
present program aims at taking care of XOR, NAND, NOR and 
NOT gates, 

7.2 Difficulty 

The gate equivalence for XOR leads to change in trGe~ 

level. Some method must be evolved to find sub-branch 

leading to change of level. Whole fault tree data must be 

modified accordingly. There are many zeros in some of 

also 

basic event entries, thesG_^must be taken care of, 

7.3 Approach 

XOR - gate has the following gate equivalence 

A(+)B = A.B + A.B (See Figure 1) 

NOT - gates are taken care of by following - 



49 


De Morgan’ s law - 

(i) a' +' B = A , B 

( 11 ) XTB = A + B 
Example 

See Figure 2 

XOR gate is first repJaced by combination of two AND 
and one OR gate. Negation at any gate is rOToved by 
changing AND to OR or OR to AND gate with input negated. 

This process continued till we come to negated basic 
events and tree containing only AND and OR gates. 

7,4 Program Details 

As earlier we denote AND gate by an even number and OR 
gate by an odd Number, XOR gate is indicated by an odd Number 
larger than 100, NAND, NOR are indicated by a negative even 
or odd integer respectively. Next, we follow the following 
procedures 

(i) First entry of each basic event input data is 
checked for XOR gate. If no XOR gate, we go to next basic 
event. Then the second entry so on, 

(ii) If any basic event has XOR gate, level of the 
branch containing Iphat basic event is found. Level of 
branch is the number of gates encountered when we traverse 
from top event to that basic event by minimal paih. If level 
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of the branch is equal to level of tree, Level (L) of the 
tree le increased by one. 

(ill) All the data right to the XOR gate is shifted 
by 1 with event number placed in array GN (I, LB + 1), where 
LB = level of tho branch concerned, 

(iv) GN (I, J) =1 and GN (I, J+1) = 2 with GN(I,J+2) 
negated, 1 indicates OR gate and 2 AND gate. 

NOTE; GN(I,J) = Gate no, for Ith basic eveit and at Jth level 
of tho branch containing it, 

( v) Next line consists of copy of GN(I,J) V j with 
GN(I, J+l) = 4 and GN (I,J+2) unnegated. Hence all the 
entries below Ith line is shifted down by one, 

(vi) Next line is looked into. If it contains the 
same XOR gate at the same level from top event, we do the 
following; 

(a) GN (I,J) =1 and GN (I, J+1) =2 

with GN (I,J+2) unnegated. 1 indicates OR gate 
and 2 and 4 AND gates, 

( b) Next line copy of (3^ (I,J) V J.GN (1+1, J) = 1 
and GN (1+1, J+1) =4 with GN (I, J+2) negated. 

(vii) This procedure is followed for all the basic 
event entry. If XOR gate different or at different level 
in tho branch concemed, we go to step (ii). 
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(viii) When aJl entries has been looked into a 
particular level from top event, we consider the gates 
at next level from top event, 

(ix) After entire input file has been modified, 
wo restructure the tree according to the highest level 
decided in considering entire input file. If level is L, 
basic event in each line is placed at GN(I, L+1) and other 
entries are filled with zeroes suitably, 

( x) Negation of a particular line is considered next. 
If any gate is negative then even is changed to odd and odd 
to even (AND to OR gate and OR to AND gate) and next non-zero 
entry is negated till we encounter the basic event number. 

This procedure is repeated for all lines in the input file, 

7,5 Output 
— 

We get a fault tree containing only AND/0ft gates 
Some of the basic events are negated and some unnegated. 

The resultant fault tree can be easily fed to minimal cut- 
set finding program to obtain minimal cutsets. These cutsets 
information can be given to fault tree quantification 
program to obtain top event unavailability and other 
parameters. Our minimal cutset finding program PCOIVCPis 
such that actual gate number is unimportant as long as an 
oven number indicates AND and an odd no indicates OR gate. 

Thus repetition of 1 and 2 as gate numbers is modified fault 
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troG doGS not cxGntG problems in further analysis. 

Modified fault tree and actual tree data can be traced 
easily from the print out. Two examples has been shown in 
the print out, 

FIN A WORDS 

We have complete package of fault tree analysis for 
nuclear systems. These systian can contain electrical system 
(o,g, scram signals). Only thing left is fault tree cons- 
truction program. But normally manual construction is done. 
Once fault tree is available, it can be analysed for top 
event failure probability using our package. Tree can have 
all kinds of gates. Final analysis depends on availability 
of basic event failure data. Assumption of constant failure 
and maintenance rates may be irritating. General time 
dependent analysis could be carried out by suitably modifyxng 
the programs* But this kind of exercise is unrealistic, 
since there is dearth of failure data. 
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APPENDIX - T 
DATA TO BE USED 


A. WASM-1400 [ i,z] report gives the following data: 


( 1 ) Pump 



Mode 

Q or 

X 

Median 

(a) 

Failure to start 

3x1 0^ 

- 3x10"^/d 

IxIO^^/d 

(b) 

Failure to run 
given start in 
normal environment 

3x1 0“^ 

- 3x10“^/hr 

3x1 0“^/hr 


(ii) Motor Operated Valve 


(a) Failure to operate 

3x10“'^ - 3x10“^/d 

1x10“^/d 

(b) Failure to remain 

3x1 - 3x10""Vd 

1x1 0“^/d 

open 
(C) X 

1x10“"^ ~ 1x10“^/hr 

3x1 0"'^/hr 

( d) Rupture X 

1x10"^ - 1x10“Vhr 

1 

1 x1 0"®/hr 


(ill) Check Valve 


( a) 

Failure to open 

3x1 0”^ 

- 3x10"'Vd 

1x10“^/d 

(b) 

Internal Leakage 

I 

O 

- 1x10"Vhr 

3x10 '^/hr 


(iv) Relief valve 


[ a) 

Failure to open 

3x10“^ 

- 3x10'"Vd 

1x10~^/d 

:b) 

Premaiture open 

3x1 0""^ 

- 3x1 0~^/hr 

1 x1 0~^/hr 
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(v) Pipe 


( 3 ) _ 3" dia rupture 

3x10 — 3x10 ^ /hr 

1 x1 0“^/hr 

(b) _ 3" dia 

3x1 0~''^ - 3x10"‘^/hr 

1 x1 O”’' ®/hr 


(VI) MTTR for 

Pump = 37 hr (typical) 

Valve - 24 hr (typical) 


B, Green and Bourne book gives the following data “ [ 3 ] 

(a) Pipe - 0,2x10"^ hr”^ [ Failure rate] 

(b) Control valve ~ 30x10 ^ hr*"^ [ Failure rate ] 

(c) Solenoid valve ~ 30x10“^ hr~^ [ Failure rate ] 

A.H, Earl has given the following information about pickenng 
A and Bruce A life time incapacity upto the end of 1981 (CAI-JDU) [ Ig ]; 



Pickering A 

Bruce A 

Heat transport pump 

0,2 (years) 

0,2 (years) 

Pressure tubes 

4,9 (years) 

0,3 (years) 

Boilers 

0.5 (years) 

2,4 (years) 

Turbines & Generators 

5,8 (years) 

6,5 (years) 

Heat exchangers 

0,9 (years) 

0,0 (years) 

Valves 

0,4 (years) 

0,0 (years) 
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J.R, Aupied and H, Procaccia have givai the following data [ ii] 
for valves s 


Sys tern 

No*of 

pper. 

nrs~ 

MTTR 

. , , ,1 

Mean j 
uaava3 * 

hr"^ 

on demand failure 

Pneumatic 

valves 

700000 

22 

53 

16x1 0“^ 

5x10“^/d 

Check valve 

600000 

39 

54 

8x1 0“^ 

0.6x1 0"^/d 

Large flow 
rate valves 

350000 

15 

70 

50x1 O"^ 

4.5x1 0“^/d 

Small flow 
rate valves 

330000 

8 

49 

1 

54x1 0“^ 

1 .7x10"‘^/d 


Steam valves 


Safety 

560000 

34 

1 ^ i 

59x10“^ ! 

2x1 0“^/d 

relief 

valve 






( low 

pressure) 






Safety 

relief 

2.25x10^ 

24 


^9x10 ^ 

8x1 0~^/d 

valve 






( hi^ 
pressure) 






Check 

valve 

1 .43x1 0^' 

11 ^ 

1 


6x1 O"^ 

0.03x10“^/d 

Motor 

1 .425x10^ 

^ 20 

i 

0,7x1 O”^ 

0.7x10“^/d 

operated 

valve 


t 

\ 




Pneuma tic 
valve 

340000 

24 


65x1 0~^ 



1x10“^/d 

1 

t { 


From A.H. Earl's data [ 1 


for pump 
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unavailability = ^ 

Q = 1 - \ = 

= 9.171x10"*^ 


8x10 


at the end of one year. 


1 


In 


1 

1-Q 


for Heat exchangers 



^ = 4.185x10' 
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Simps 


Error 

S 

Hence 

e = 0, 


APPENDIX II 

algorithm of numerical integration ^^ 3 

on’s Rule f 



" i=i! 3,5, .,n-2 I ( 

= § ( f,+4f2+2f3H-4V2f5+.. . ) 

f IS to be tabulated at odd number of points. 

If accuracy upto Nth decimal place is desired 
5 * 10“^. 


(11. 1) 

(11. 2) 
(II .3) 

(II.4) 
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APPENDIX III 
ITERATION TECHNIQUE 

Bisection Method [ >'5] 



We begin by picking two trial points which enclose the root 
this IS indicated by being of opposite sign. 

The interval (x^, ) is bisected and mid point denoted by 

X 2 » i . e . ^2 = (^1 + ^ ^2 

If f(x 2 ) > 0 then root is between and X 2 . Hence replace 

Xjj by X 2 and search for root in this half interval. 

If f(x 2 ) < 0 then root is between X 2 and x^ . Hence 

replace x^ by X 2 and again bisect the interval. 

This bisection procedure is repeated till search 

interval is smaller than the precision with which answer is 

wanted. 

Note that this raeti:^od always encloses the root in the 
search interval and search interval is halved each time. Thus 
in 10 iterations, the search interval reduces by2^° ~ 1000 

and in 20 by 10^® ~ 
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2 

? 

2 

3 

3 

3 

3 


4 

4 

4 

5 
5 
5 
5 
5 
S 

s 

5 

7 

7 

7 

7 


II 


23 

24 

25 
22 

23 

24 

25 
22 

23 

24 

25 


4 22 


23 

24 

25 
?2 

23 

24 

25 
22 
23 


24 

25 


I 22 

8 23 
8 24 
8 25 


22 

23 

24 

25 


22 

23 

10 24 

to 25 

Ij II 

U 2 | 

2 5 


12 23 
12 
12 
U 

13 


! 2.56 

at' Sets 


7 


u 

24 

u 

25 

1 4 

22 

u 

23 

14 

24 

14 

25 

J 5 

?2 

» 5 

23 

15 

24 

15 

25 

IS 

22 

1 5 

23 

IS 

24 

16 

25 

U 

12 

17 

23 

W 

24 

17 

25 

IB 

22 

IB 

23 

13 

24 

18 

25 

19 

22 

19 

23 

If 

24 

19 

25 

?0 

20 

II 

20 

24 

20 

25 

ll 

II 

21 

21 

II 


,fCJ38/8Ar:H 


[ 15100,15006411 

Lo^isd-oCC ri-irU6 at 14:10:55 do l3-Dat-i35 
Raotiga; 0:00:0?, KCS:45j cpooect tine: 0:03:07 
ols< R8a5«: 235, 'frites:li 




’ Exonmyle 1 

^»concp 


X 


BAirCON V3irsloa_lt>2(20 


input* £r:).Ti dskc 
DUCpUt "" 




ss;ciu3hc:2l 2752' in strea.n I 


„ tD; OSKC 
jipof paraiietelrs 
rt'neirOOJ 31100 Coreiioop 


UAiqu.eiYES R’sis’tart'iYES Output*: 


'tOG-IlM J 5 i 00 / 


jorsi VI V] 

J 3 P. 9 tn:eir 

1340 29.^10 


^DEFER/5P3OLlAljCi/iTlMS!:^50/COREM00P/tOC:A;CE;:i0/?<AME4'':M4NaJ KOMI 
ivANPUK OV 3 A(vj i ill 15 
3 : 0 bs same PPM: 27 ] 


[f5„r6fg?s 


t 020 C 


i 


85IS8 


80 ( 

mi 

JMS 

mn 

01700 

01800 


8 


mn 

«li 

111 

0290? 


»480( 


V "86 


sat* 


PROGRAM MAME: D^OOK'IJT?' 


THIS PROGRAM OSES KITT FORWATlON FOR FIMOIMG SI^rE'Ji UMAVAI-' 
LIBILIiy AMD SHORT CUT ClASiClOLiArlDN METHOD GIVEM- Biit j.. FUSSStt 
FOR FIMOIMG.. SISTER PAgAMEffERS LJKE Wpi AMO' b'AMiBO'AC SJ , 3, tO 


AND V(Oit') ARSI ACCURATE^ FOUND' FOR BASIC EVENTB.r 

TKf B<T ur»» Mr«. c'Ve 'I's'M nra-^ ais.ii 'ifi/PinT Air* tSftrtriB'rvun t!* li alIa 


IM fInBIMG. STSTEM‘'Qr33 3LRA.C^.SirrMG''PRDCE0UR'B”HAS BEEM OSgD.* 


FIRST. BRACKSr gives D(S)MjIN' ANO' SECOND Q(S’)WAX. 

FOR BASIC! events Wfi}^ SAM SP.SCli:Fy M OUT Or M S^STEHlt'iN THE PRO-: 
‘ ■ ■ "* :liarli me: can 

- -C. SATES rfE 

_ _ _ _ OUT OF M 

GATE 

t 5ERE ARE THREE OPTlDM- TO THE: USER.FlRSr IS WHgM ABB' CUTSETS 
ARE TQ BE! TAKiSH INTO AC;C1,6'SCDM0- IS WHEN..MAX*' ORDER, OF CUrSETS 
INTO ACCl, THIRD' IS WHEN MAX.' PERC'SNTASE OF THE 



are to be: takej 


ni'K* *w WD- 1 , A w ’£ uji-rvv' 

MAX OF AU* CUTSETuQ' IS SPECilFISD,' 
TOD “ ■ ■■ “■ 


-1 METHOD, OF, l:M.PUrTlN,G THiSBEi I:MF0RMAT;I0N: -1 
IN FitE.NIAME FQR2.4 .OAT G.IVS1 


IN THE FIRST BtlMSi TIME AMO' DECIiMAL ACC:URAC:y OF WtD.t) AND' VCO,tO: 

wanted, 


IN SECOND’ BINSi C3MP. MO, UftM’310Ai,MfiUi.M,N,MN 
M AND N are: as MENriDHED' AiBDVE'b'WHBN ABL 


FOR ABC.! BASIC. EVENTS. ' 


the; BASiCi: evemts are 


OVER LAST' BINS! MUST BE EIIITHERI A1 BBANCK BiME; OR A' LUME; SrARTIMG 
W It HI 0 * 

MN=0 PlsBAMBDAl ;P2aME!J; (O'AT'A; TO BE GIVEN) 

MN NpT=0 => Plx2|r) ..pIls^ID (DATA TO BE Gi;VSN)^. 

NEXT line; HAS QPrlOMsT 3P.r’lDN«0 s>TAKE alb CUTSE'TS.' 

QPTIQNsji =>TAp «AX..3BD|R! ,SPEC!l,raATi:OM. 

OPTIOHsl OR' MORE 5>MAX | IS' ASSIGNED, ' 

IF opTiOMsi MSftcT LINE: c;5MrA!r:NS‘ MAX orders? 

IF OpriON»2 OR MORE! IT CDNTAHNS ta? 

NEXT LINES HAiVfe MpTMAL' qU:T®ET INFORMWIOM,' . . _ • 

these: again must end MlT'Sl SI:TH£;R: a BLANCK: lilNE: or a U<INE 
STARTING WITH 0 

LAST line: HAS M ANDi N’ of: TOP* GATE;. 










W- ‘in «ii iii'w . 

LL*' LOWER* limit OF 
UL- UPPERVSlMJT^Of 
CN- component: NO.V, 
LAMBM-: mean L'AMBO'A- 
ttfM- MEAN W 


UMAVA!ILtABEl.lTy 

D: : - 




3 


86000 

85100 


lun 


HIM 

97500 


97800 
97900 
98000 
98100 
) 8 ; 


>37( 

9890( 


SUSS 

ISfSS 

I8ISS 


10500 

ifiSOO 

IO 7 OO 


11 


00 



>si 




20 

25 

30 


35 


35 




cw(ioo5- caTss:r « 

KSCIOO)- curs SI ORO€R ARW 

CfAMBClOO)-! CUtSEl DAMBOft; 

SVC 100)- vet) 

F- Q By MARKOV'S FORMUBATIUV 
£1" 

F2— V ( t ) t m 

ISDICAISD' EARLIER 
P2- ME)i 3R wer') ,, 

T- XIME 

Ml:.°i"ouT^bFi"|ysra 

k^- INDICAXOR WHElaSR ?!■ IS' LA«BOA OR Q? P2 IS ^SU OR « 
DSPEMOUSS. 3S WHSirUSR' IT IS 0 Or NOMZSRO 

ft'' ft 


• W ,«# .*# liti!** m M 


REAL LL^LA«BDAL MSP, LAMS'* 

iMTEGER :?i.cs»^pri5u,C3Jvr,PAC!r , 

DIMEMSION 5 C 100 ),SwCiO 0 )#:B(lO 0 aOO),Ca(lOO),C( 4 :ClVa).Ki>CiO 0 ) 
1 CLAM 8 ( 100 ),SV( 100 ) 

external F, FI, F 2 ,FACr 


COMMON P 1 ,P 2 

F§RMAT(f 8 ,l!ixll 2 ) 

WRIXEC 5 # 25 ) 

FORMAT(/,lX,':OMp. fO t' . wCt) 

QC?r* 5 '' 

DO 100 1 = 1,100 ^ ■ 

REAO( 24 f 305 cNiP 1 ,P 2 , 4 |,M, MlS 
P g R M A T ( 1 X # 1 1 4 ll 1 ^ 1 si 0 • 3 ) , . 3 '( 1 K , 1 2 ) ) 

lF(CN,EQ.<hGDrD 1 05 
lFfMN,N 

swr 

S V c 

eiS:tirrr4:S::;l:!!i:8 

ESig’ir’' 

QCCN)sPl 


V'C tU 


SW(CN) /(l.^aCCN)) 


||SigS«5 

PlsLAMBOA 
P2s0, 

CALL ET<F,3CCN),N1) 
MEU-P2 . 

!?fS?^5rS5gJt§^5S’ 
pclJi i 

li2»N 

DO 40 J»M>N 



10 


•ST)»*CM2-j;)) 


Jf) ♦♦(«-!) I ♦SVC CN) 

if , 

Ept|||swiC!l3/Xi ,'-:2tCN3> ' '■■ 

ife' 


ORflOM 




13S00 
1 37 DO 
13300 
13300 
14000 
14100 
14200 

11122 

11^22 

llisi 

14700 

14300 

14900 

tSOOO 

15100 

11122 

llin 

15500 
15500 
15700 


120 


j |300 


300 


ffio 


000 


20 


filoo 

mil 

mil 

Boo8 

17100 


iim 


► 0 
17400 


im 

nu 

l? 3 S 


8400 

8500 

8SOO 

8700 

3800 

8300 

9000 

9100 

l?22 

9300 


‘l|g» 


1 


tflll 

9900 

jgjgg 

¥ 

ii! 


130 

HO 

145 

143 

145 


ISO 

ISO 


1^5 


170 


220 

230 

250 


250 : 


2 S 

21 


25 


2S2' 

|B5 

w\ 


290 

01 


34 0 


FQRMATC10X,r2) 

G 0 io’l 43 

fogSi|t5itgJ!r, 

MOsloO 

«SJ?E15».146) 


DO 165 1=1,100 


t 




DaiDia+i t) 


0 »(t) ' 


READC24» 160) 9 (CS ( 1 , 1 1) , Ils.i,, l OQ) 

C> n r»l 111 Ji m ,ir i »% A # ,J, tV 3 n C *» ' m r ^ m ^ ^F 


_ i f A o y i , . , , , 

imiH 


17 0 


COUNlsCOUNr+I 

continue 


iCUlSET lambda E:aOAriDNS.' 


KSCI)»0 

;) = i, 


Hioo 

).E3.0)G3r3 205; 


CQCI 

C«CI)=0. 

DO 200 f 1= 

kinUiiU. 

CQ(I)=!5Q(I)^2Cj) 

DS 520 Lsl,KSCl) 
IE(0.|Q^C:SCI,Li)) GOTO 220 
KsCStI 


CONTtNul 


C) 


C»(I)sC«CI)+S4CJ)»S 

continue 

CLAMB C IlsCM (I ) / C 1 . -iCQC I) ) 

CONTINUE 

AS PER OPTION ACITION IS TAIKSlN.f IF OPTlONsO , ALL CUTSETS TAKEN 

GREATER' THAN MAX ORDER 

SPECIFIED- THEN 3 IS REDUCED* TO ZERO, 

IF 0PTI0>|>=.2 IF ACCORAC:! CM]AX 3-3) /MaX Q *100 >% ASSIGNED THEM 


AGAJM_Q_R£OOCf 


ocsb^ra^ZERo.f 

S. 0)3313 


SS“Sto 4x2, COUNT 

I £ C C 2 U ) , *gI .Cl'^ ) CiVsClQ ( I ) 
CONTI 


)»l50./cr) .LEI, PC!) GO TO 270 


CONTINUE 

DO 270 Ixl.COpT ^ 

IF(OPTlON,fi,Q.l )GOT3 .262 

GOTO 264 . 

IFCKSCT),LE:.M3)G3T3. 270 
CQCI)® 0 t^ 

CLAMBCiJsO.t 

CW(I)eO, 

CONTINUE 

WRITE ( 5 , 58 2 3 , 1 1 # r» Cl'^ ( I ) 9 CSiAiMBC i;) 9 CQ 1 1 ) ) 

FORMATCUf l2,lX9F8.4 9 3(U,sa3..6)) 

ilpLATlON OF 5YSTS?< 0 MAlX,a' MllN,0 ME AN 9 LAMBDA, w( n 

LOOP’fiNdS first TERM] OF SYSTEM 3 EOUATION,- 
DO 290 Ixl. COUNT 


88 Itg'SlP-**' 


SECOND TfeRM. OF SYSTEMi ,.3:-SaUAXION AS 

cour* 

rl-l- 


:dri ,» 



i;{) 


I f ceiii £|?!2;5 m mp\ 



21300 

3 '2 0 

21400 

21500 

3'3 0 


34 3 

ims 

V i 

il?SS 

Illoo 

lil|S 

22500 

22700 

350 

22800 

22900 

23000 

360 

23100 

lliu 

23400 

37 0 

23500 

23500 

23700 

375 


continue 

cori=cori.+:um 

continue 

CONTINUE 

LLsULtCORI 

cgR2=6. 

D0^400^1”3°''3JmT STsrSM 3' SaOftflON AS C3R2< 

DO 390 J = 2M-1 
DO 380 K=1,J-1 
CUMsCO(I) 

DO 360 Kl=ifKS(J) 

Dp 350 K2 = 1,((.S(I) 

C^TINU^^^ ^ ^ ^ ^ ^ 

KlisCSCJ.Kl) 

CUM=CUM»3C KL) 

continue 

Qp 378 Ko=1»K;S(IC) 

00 370 K:2=l,iC:S(l) 

IFCCSC K, SO ) , S a.C S C I.r K2 ) ) SOTO 37 3 

continue 

DO 375 Kl=l,SS(J 

“ "oj.eo.c 


IFCCSCK.SO: 
CONTIMUe 


(Ji, Sl))30r3 37 3 


2||23 

23900 

24000 

ills 

mss 

mss 

illoo 

mu 

25200 

25300 

25400 

2S500 

25700 

25805 

mu 

ills? 

illSS 

UlU 

26500 

26700 

26300 

2/ 000 
27100 
27200 

IS 8 

ifiss 

27700 
27800 
2790a 
280 
Ifi^O 

Illoo 

IllSS 

tooo: 

fOOO 


37 8 

36 0 
300 
400 


450 


460 

465' 

Cl 


470 


475 

480 

490 


n\ 


'’1 

j; 


KpsCS(KrKO) 

cgM=cuM»ac sd 

continue 

COR2=COR2+Cy« 
gONTlNUE 
continue 
continue 

SMisfeRr?u|#LD<) 

READ(24f 450)M>N 
F0RMATC2I31 
IFCM.EQ. 0)5013 465 
AS PER OUR FORMULiAriON r 
STsQM 
QM=0, 
srisUL 
ST2=LL 
UL= 0 . 

IiLsO 

DO 460 J=M,N 

I vir ACT c fi 2 > / c FAc:r( J 2 ) ♦f 4 cir c n 2 -ij :2 ) ) 


TO. FIND' Q FOR M' OUT OT. sySfEM 



ONTINUE 
NMsO. 

LAMBMsO. 

LAMBM=EAM8,i+CLAM3i(i;) 
5 ?? 5 ii§fo) 50 r 0 : 475 


RT=F LO AI ( F AC I C .N ) /% f ACT ( M -jl ■) *l ACI J » -M ) ) ) 

W M=RT» C STt * C M fl 5) » NM /1 1 • <- C ST 


LAMBM® WM?i( 1 #*3 m 
NRITEC 5 f 480 ) . 
FORMATC/ZI/IX, ' 

1 QCS)M 


>wm) 


LAiMBDAlS') 


CS)HEA;!I ') 

f E C S J 'I^O ) C r, L Af B« , WM , piLf. m, 3M ) 
format C 1 X f F8 .1 , 5 C 1 X , El 3 •6) 5 
STOP 
END „ 

||•EKcA[^^s|igMifaS 

OUTPOtI Also SfER. USED FOR! OfiBigEO ACCHRACX „ 
FUwcTIuN^ro 8£! fed THgDOSai Af EXTERN a|# FONCIION 
ri DEFiilj; S The: EtiM?a»roN: *' 


»Csj 


Q(S)MIN 



29303 

29503 

29600 

29700 

29800 

30000 


RC.AL lOfll 
COJ^MON &1/P2 
H=CX2-Xl)/2. 

1 = 2 

Sl = Fl(Pl,P2,Xl)+rl(Pi ,P2VX2') 

S/asO - 

S4=Fi(Pl,P2,Xl+H) 


3«M 

loioo 

30503 

iini 

I 8 i 3 § 

31000 

fllOO 

31200 

31300 

31500 

31500 

31700 

31800 

3I9O3 

32000 

32100 

32200 

32300 

32400 

32S00 

iiiSS 

ims 

um 

im. 

iiip 

33900 

34000 

34100 

34200 

34300 

34500 

34500 

34700 

34300 


li=C§l+ 


!.!S4)nH/3.0 


36100 

l§29^ 

35300 

35400 

35500 

35500 

35700 

35800 

35900 

11888 


10 


5^ 

5' 


tWf I 


999 

1000 




1050 

1300 

1400 


lF€Il.Ep;o;-)S3f0 6 
I F C A B S ( C 1 1 -;l 0 ) / 1 1 ) . L S , 0 . 5 1= 1 0 . * * ( -:PI ) ) G 0 r 0 5 

54 *52 4" 

S4 = 0. 

X=Xl+CH/2.) 

DO 2 Jsl,I 
S4=S4+F1 (PI,P2,X) 

X=XtH 

COfjriNUE 

H=B/2. 

I=2»I 

10 = ll 

11 = CSl+2..*s2 + 4.»S4)»(H/3i.') 

GOIO 10 

CONTINUE 

RETURN 

END : 

¥!!!¥$**$$$*$*$$$***:* ¥*^**^*******4: 
FUNCTION FCP1,P2,T) ^ ' 

IF(C(Pl+P2)*r5.CT.50.-)F=Pl/l(?ltP2) 

I F ( ( C P 1 + P2 ) * T ) . U e:. 60 . 0 F =.( PI /I PI + P2 ) ) =•= C 1 . - E XP ( - 1 * C PI + ?2 1 ) ) 

RETURN 

END 

$*¥*¥***********¥***¥¥¥*¥¥¥¥¥*¥***$¥¥*****¥*¥*¥¥¥* ¥*¥*¥*¥**** 
FUNCTION Fl(Pl.P2#r) 

IF ( ( ( PI +I2 ) * f ) ■ ul;# loll fi =1 c P 1 ipl ) Ic p 1 ^ *2 ) ♦EXP c -:?♦»: PH-P2) ) ) /:( pi +p; 

RETURN 

t¥¥¥‘¥***¥*****¥'^¥*¥*¥** *¥¥¥¥¥¥¥***¥***¥******* ***''^********¥*****^ 


If ( ( ( P I <-P2 3 ♦ r ) . L El. 50.0 F 2‘sl ?A ♦ P2 / C Pi +P2 ) ) ♦ C 1 . - 


EXPC-C*tPl+P2))) 

REfURN" 

END 

^*¥*)¥i¥¥:*¥*¥******* *****¥¥*¥$¥*¥***¥***¥*****¥*****¥******* ¥*¥*¥* 

integer FUNCriDN FACTCN) 

FACT—'l 

irCN.EQ.OlGOrO 1000 
DO 9§9 ial,N 
fact=fac:t^i, 

CONTINUE 

fill* .1 I * t I’ r t- 1 * I t ». ♦ I f f i- r-»4 t ,1 I,., I 

SUBROUTINE EfCF,S,N0 
COMMON P1,P2 
11 = 1, 

YOsO,- 

if(ASs(cn-y00/yi0,Gs;.'0,5#i0.<^*c-N3)GOxo 1500 

'■::P 2 Sp'X ]f;t + I ■ V 'J /'iZ » : 

Z=FCPl,P2fr3 

lF(O-F(Pl#P2,r)3l3OO,153O^l400 

?lap2 

goto 1700 

l0eP2 



15»3 


36&1>!) 

ill:?????]®” 

e-io 

2 i;303E-l2 

3 fl.30!)E-lO 

t :};s^o”3l:Ji 

n-Mipi 

l| liUlUi 

8.3.00E^:0& 
3,50O£;-O6 

2,090S-.39 

mmu 


GOTO 1050 
RETURN 

end 


ii 

13 


I 

17 

18 


+S,000E- 

4.000E- 

4-5.000E- 

+4.000E- 

+1.000E- 

+S,OOOE- 

tl.OOOE- 

+5,000E- 

i-.mt 

Z.OOOE- 

5.000E- 

4.000E- 

5.000E 

4,000E 

2.500E 

|•000E 

IrOOOE 
5.000E 
5.000E 
5.000E 
1.200E 
l.OOOE 
2. 500E 


“I 

n 

04 

0 ? 

ii 

09 

§i 

U 

-SI 

"04 


il 

■05 

- 0 | 

-07 

-08 

-04 


3 

2 


2 

2 


4 

3 


3 

4 

3 

3 


I 

5 24 


6^ 

7 

7 

I 


2? 

I 

it 




.i 

I 


:■ 11 23 

11 24 

11 25 

12 22 

12 23 

II II 
II H 
li II 

A li 

11 II 
, It II 

: 15 24 

15 25 

It II 
It II 
II II 
II II 

I il 

18 24 

Jl It 

l| 22 

19 24 

II II 
IS II 

29 25 

21 22 

i II 

21 25 

3 

3 3 

,e:x 

h 3 S 4 i/l 3 
LCMK.: ^ Loading 
[D'IKx.i/ 4. S’xecutlonl 



9 


04iP 

f 0 t 

wCt) 

V ( t ) 

j|tO,t) 

** 

350103.0 

0,999965E-:1D 

0.175344E-:i2f 

3.(350394E-04 

7 

350400.0 

0,100000E-11 

3.140216E-:15- 

3.350400E-06 

3' 

350400,0 

0,999965E-!10 

3,175344E-12' 

0. 35 03 9 4 E- 04 

4 

35O490.,O 

O.lOOOOOE-11 

0.140216E-:15- 

0,350400E-06 

5 

350400..0 

0,787257e.]05 

3.30751 5E-;33< 

0,003000E+00 

6 

350400.0 

0,999965E-10 

0,1 750 44 E-1 2 ; 

3:.350394E-04 

7 

350400,0 

O.lOOOOOE-ll 

0:,1402i5E-;i5' 

0 , .35040 0E-:06 

8 

350400.0 

0,999965E-;10 

0,17S044E-;l-2i 

0.a50394E-04 

9' 

350400.0 

O.lOOOOOE-ll 

0.1402lSE-:l5' 

0.350400E-06 

10 

350400.0 

0,681855E-!07 

0.83OS13E-:35i 

O 4 OOOOOOE+OO 

11 ' 

3a04»0,0 

0,359U8E-:05 

0,171662E-0.4' 

O.OOOOOOE+OO 

12 

350400,0 

0.999965E*!10 

0.175044E-12! 

0...350394E-04 

12 - 

350400,0 

O.lOOOOOE-ll 

0..140216E-:i6' 

■0,350400E-06 


V(0,ti 

0.'3D6769ai-07 

0.245 5S7ai^ll 

0.3 36 7 69 £-07 

0,245SQ7£-11 

0.'000003£+03 

0.306769£^Q7 

0.'2455S7£-ll' 

0,305769£-07 

0,245Si7f-lt 

0,000003 £+00 

0.030000f+00 

0.336769£-0/ 

0,a45&67£-U, 


[jA'lBOA 

O.IIOOO03E-139 
O.lOO30OE-;ll 
O.iOOOOOE-iOJ 
0.i00003E-;ll 
0.8 04236E-;3S 
0.4 30003E-i39 
O.i 30a03E-!ll 
O.'1303OOE-;39 
0.l000 03E-;ll 
Q.633278 E«:37 
O.'39O094E-;35 
O.lOOOOOE-lOf 
O.I 3 OO 8 OE-II 


HEU 

0,500 30 2E-08 
0.400i6lE-l0 
0.500302E-03 
0,400164E-10 
0.i45657E-Ol 
0,500002E-08 
0,400t64E-ia 
0,500002E-08 
0.400l64e-l0 
J>...'398929E-12 
0,216182E-03 
0,500002E-08 

0 , 400 l 64 r -10 


at t 

0.330 38 7E-; 
0.350397E-i 
0,350 387E-: 
0.3 50 39 7 E-: 
0,211l.2 3£-: 
Q,a5038 7E*; 
0.3503'97E- 
0,350 38 7B-, 
0,330397E-: 
0.238 211S- 
0,9943S5£-; 
0,3503S7E-; 
0,a50397E-; 



14 asOtOD-O 

0.999965E-:10 

5.175044E-121 
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ir-i 15133/ 
31 I I r 
ijs?. otn.eir 
? l3-:De 

0,3 i/e soT'e 

CT.F3R 
S38S23J 
03 :: 

03 C! 

03 Cl 
03 

03 :: 

33. 

03 

33 

03, 

.03 

03 100 

33 
00 

■03 I'S 
i03 200 

^03 
'03 

too 

103 20 

133 
I '03 

203 25 

JOO ■ 

403 
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503 500' 

700 

303 30 
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183 
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803 
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PROGRAM MAME MlMCELS 

PROGRAM FOR MIMUfftL- COl AiM3' P.ftirH SETS FOR ft SLSCrKICftL 
integer CD 

common CDC7»53) elMMC50) , IDDJMTCSD) 

C0MM0N/PftTH>MRrh5 , 43 ) J 'm?,T I c 25 , 403 , MpTO C 25, 40 ) , ,M?.U 0 3 , 
OIMENSIO^I rS03E(l0) 

READC22£l033M30E 

F0RMATCI3) 

00 20 IRsi,7 
iSsJ 

READ{22,2O0)C:dCi:R,IC),IC!=-i:5,i:s<-l9) 

F0RMAT(23I3) , « „ ^ „ 

I F ( C 0 ( I R , I S + 1 9 ) , EiQ . 0 ) GO r 3 2 0 

ISsIS+20 

GOTO 15 

continue 

DO 25 1=1,7 
DO 25 L=l,50 

IF ( C 0 C I , U0 ,€3 .0 3 CD ( 1 , 00 = -11 
REAO(22»l0g3MDS „ 

READC22, 2033 (INODE CD) , 1=1 ,4 OS 3 
WRlTE(5,30D)a:N03€(D),Lsi#M3S) , , 

FORMAT(/ , lx, ^Source; ,»I0DS5S' ' ,23133 

DO 30 1=1, MOS 

INNCINODEII:))*! 

CADD TRE£;(MOD2!3 

call Bcur 

STOP 

END 

&&&&&&& Si 

SUBRQUTIME TRSECND3€) 

|NTEgER.C,9. 

,43),,M?,11O0, 

LEVaO 

MPT(l£l3=NODE. 

DEyaDEV+l 

IFCMPTCDEV,1).EQ.03GDT3 73 
Dli®0 
DNaO 
LLaDD+1 

IF(MPTCDSV,LD).E3. 3)3013 20; 



"MAMOJ IC D M A R " 

SC stem 

15),MP,DP 


'*S&&iSifi&SiS«Si 5 «& 
1 5 3, HP, Dp 
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i)? 

'J*) 

0 0 
53 
-35 
55 
135 
• 55 
>33 
13) 
>3) 
too 
100 
>35 
>03 
100 
333 
900 
000 
100 
203 
305 
403 
503 
603 
703 
300 
903 
300 
103 
1203 
1303. 
1403 
1505 
(503 
(7 03 
5300 
5903 
?0O3 
>103 
>20 3 
)300 
9400 
9500 
9800 

im 

9900 

0300 

0100. 

0200 

0300 

0400 

0500 

0500 

0700 

.0300 


32'^ 

35 


40 


45 

43 

44 

45 
47 


13 

50 

52 


55 


7 0 

1000 


33 
35 
90 
10 3 


1010 

U5:^ 


lFCIf^NCH?TCLiSV^»LiLi)5 .£0.4 >3313 33 

Lf4C = p 

MPP=MPjCLi€V#tt4i') 

lF?IMl?C3CL|::.!4PP)).f24:>^|P 

tF(CD( LN^Zl. ■’<PP) .1^£. -.1 )'j34 3 *3 

‘4PTDCLe(r,LU)=!^N 
GOTO 30 

IF(LSV.E3.1)33t 3 Sa 
iTTsl 

LLI = ^4PTUCLEf i 

IFCCQCLNZ# ^PP3 .SO* '4PT(l<S/-llJr+t » ir))G3rJ 35 
GOIO 52 

lF(MPTDCLC\f-irT#Lijr“lM) .GJ.OOSOTO 48 
iHsIM+i 

IF(14£.T.E2.IM)30T3 43 

DO^io^lTsMPlO ( GE4f-I.Xr . L.L|r-||l) n » n"35^ * 

IFCCOCLI^C f ^PP) ,EO.'4 PX‘l»*''^"!I 4'T+1»XT) )G3T0 35 
iTIsITT+l 


iF(iTTlE3,LSV)G0I? 55^^ 


M P P T s M P T £ L £ V -I .T T f tj LX ) 

ip| C LEVM I L|5 ) =C3 C , >4P ? ) 

MPTlKLE'^+'l f 

3 :. 4 iSl::! «NI.«AL' PAIH: FR 3 M- TRE£..;-<-W..w:.:-.r- 

.O0£r'>13,//*' basic P4T.S 4 '5 

? 4 PsO 

bo 90 I = L€V-l 

N 5 M + l„ 

GOTO 80 ^ 

continue 

RETURN 

NPsNP+l 

IjPsI 

hh»l 
LNsN 


IPalll 


20145 


=LSI“Hpl5ct,t,U.) 

LL»LL-1 

LPsLP+K 

GOTO 110 
END 
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£>£> 
93 
33 
33 
^33 
.33 
'33 
?33 
?3? 
?33 
13 3 
233 
33 3 
433 
533 
533 
733 
333 
933 
333 
133 
233 
333 
433 
533 
533 
7 03 
'833 
i933 
‘303 
fcl'OO 
t2’33 
U03 
1403 
1503 
1503 
4700 
1803 
1903 

mi 

52 09 
5300 

112 ? 
5503 
5700 
5300 
5900 
6000 ■ 
6100 

I1E§ 

5400 

6500 

,6500 

.6700 

.6300 

16900 

1|?S? 

m 


SUBROLITI^E B:J:r 

integer :o 

C0M«0N/P^^^/^?TC25,4^),«?J’J.C25,43) ,«prDC25,40) ,H?.U03,151 ,NP, L? 

cg,vimo(^ goC7,533 ,i;n5c5o 5, rroo.Nnso) 

DI"!ENSI0N LFI8SIC2 0) ,LSS::iC l3,2) j,LXrflRDC50,3),L?00f<t83.|) 
DlMEMSl3.i Lr(l5),I.NPC30/5), rNCCaO) ,MtC 25) rLMC45) ,&2ril3,20) 
OI^.ENSID'J ^31.(25,43) 

FI.RSr ORDER OUISEI'S 

M 1=0 
^iT = 0 

DO 20 1=1, bP 

lF(ICOtJNrCMP(NP,i:)),E3.'I?0 3,0r3 13 
HT=NT+1 

INC(?4P(N?,I))=-Mr 
bT(NT)=M?(?Jp, D 
GOTO 20 

lNCCMP(Np,I))=-23 

Nl=NI+l 

LFlRST(NU=MPtMP,I) 

CONTINUE 

FORMlT!f|?'^FS.sf^3RD€^^:0ls|:r i ',2013) 

DO 50 1=1, NP 
NI = 0 
M0=0 
N0=N0+1 

lF(MPCl,JiD),S0,0)G3X3 4 3 
lFCINC(MpCi:,N35).i'r.O)G3n 44 
NlsNI+i 

IfJC(MP(I,M35)-l 
MP(I,NO)=0 
GOTO 40 

I F C I NC C p C U M 3 ) ) ,f 3 , -2 0 ) GOTO 4 5 

I N p c I , IN ::( 4P ( I N 3 ) ) ♦ -a j h 

HPCI,NO)sO 

GOTO 40 

DO 50 K=1,NI 

MP(I,K)=MfCK) 

CONTINUE 

NM=0 

DO 55 jsUNf 
NMsHM+l , 

LMCMM)=l4rCj3 

DO 60 J=l,50 

iFdNCCJolNE.DGOTO 50 

NMsNH+1 

LM(MM)=J 

continue 

lC£l£i--:Ki^-:EEslc:3ND^ OR'OeiRl 
N2s0 

00^80 1=1, NT 
NL2sO 

DO 80 jsi + ltNN: 

DO 75 KsifMP 
M=0 

IFCINP(K,I),E3. 1)5313.7 5! 
iFCINCCONCjjl.GT. -0)3313 7 3 
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1-7 430 
17533 

79 

17 503 
17703 


17003 

17 903 

75 

1-8303 
18! 09 


1-8299 

18 3-0 9 


18433 

18599 

33 

1 8503 


18709 

18303 

201 

18909 

19909 

90 

19109 

202 

19209 

95 

19 309 
1-9403 

1 

j 

19500 

19509 


19709 


1-9303 

19909 


20909 

20109 

10 0 

20209 
20 309 
20 409 


20509 

20 509 

10 1 

20709 

99 

20309 ' 
20909 


21909 

21109 

102 

21209 

21309 

21409 


21509 

21600 

103 

21799 

21309 

104 

21909 

22909 


22209 

105 

22 309 
22409 


22509 

22509 


22709, 


22800 

22909 

;:'iio 

23000 

23109 


23200 

23309 

, 20' 3 

23400 

120 

23500 

204 

23500 

130 


IF C I NP € K , I M : C C J ) ) -:1 ) . '23 , <l ) 30 r 3 7 5 

H •• M 4 1 

IF(HPCK, ,eD,0)3Or0 83 
|gCgPf^^).e3vLM(J))33r3 75; 

COfiXINUE 
N2SN2+1 

L§ECtE«'2,l)=LrCl) 

IF(N2s £ 3.0)33X3 95 
wRITE(5f 201) 

Fgft^ATcijCj. '^sssDjjo 3RDE8 :!j:rssr ; •) 

write (5 J 20 iy L LI, C^SECiC LC. , I ) , uSSOX LL , 2 ) 

FORMAT(l4XfI3,‘; »,2I4) 

IFCNP, LI, 3)30X3 UO 

«;j.g I ;R 3, JRrjgR- OaiS'SI-;— -• 

n 3 s 0 
DO no i=nfir 

^IL3sO 

DO 110 Jsl^l ^ 

DO no Ksj+n«M 

I F c I C 03 M r c hH r;) ) f i :pus i ( l« c jd ) f 1:0 un i c lm x k 3 ) , lx, ) got 0 113 

M=0 

NsN+1 

IFa 2 I(nff).E 3 . 0 )S 3 T 3 Ul- 

lFa2ia,!^). £3 , L ^ c J3 . 0 R . C ,8 T :c n ?l) . e 3 . L M C K ) ) 30 TO no 
GOTO lot) 

IFCJ.GT.NDGDTO 102 
NMTsO 
NfMTsNMT+l 

If: {;2iUrp«J.Eg.9)S3.1 
GOTO 99 ' 


J3T3 102! 

I F C t, 2 1 ( J ; T ) . E 0 ; X< viX <5 ) 30 TO 1 1 0 
DO'lOS'LsI.SP 


i/ w ii V w * f r _ 

iFc iNp ( L,iy , £ 3 . 1 ) Goto 1051 
IF? J.GT.MTlGOro l 03 

lEC iNPCL , iMC { c jy y y .sa.'i ) soro 105 


IF(K.GT,MT)G3X0,104 
iFcr " 


:iNPCL,lMCCLMCiC;))»-3l).€3‘.il)30r3 105 
MS,M+l 

( I ) ? 3 si i ?'C ulf '1 ) . £3 . i* M C F ) ) GO TO 1 0 5 

GOTO 104 
COSTINOE 

N*3sN3 + 1 , , , 

GTHIRDC 

LTHIRDiN3f2)aWCJl5 

LTHIRDCN3r3)=3MU:5 

L3i(l|NL3)sL«Cj)n0OtUW (,<3 

FORMATXiXf' T' 41 RD' ORDER OaTSEir ') 

WRlTEtSf 2045 loi. C La'fl|.RDC SIiO , Kliel #3) 
FORMATCl4X,I3 »*: n 3I4) 

RETURN 



5 


2370& 

,'Ty rDR22!,OM 
tl4l38s31 J 
t 


a 

3 


i 

5 


1 2 
5- > 
7 9 
13 


2 

A 

6 

8 


END 


5 3 5 4 479 

7 6 7 10 3 8 11 

0 8 10 12 3 10 

0 11 11 01112 


1 

12 

,sx :,r,foR 

h4:58!3l J 

Li0 3'aing 

[L.^KXCr Oil’ eixecutionJ 
SOORCE NODES; 12 
FOR SINK nods; 1 
THS; BASIOi MINIMAL PATHS 

basic path 

BASIC PATH 
3AS1C PATH 
SASIC PATH 
SASIC PATH 
SASIC PATH 

basic path 

BASIC PATH 
BASIC PATH 

Basic path 
basic PATH 
BASIC PATH 
BASIC PATH 

Basic path 
basic path 

FIRST ORDER CUTSET ; 
SECOND ORDER CUTSET ; 

1 : 11 


THIRD ORDER 


2 ' 

::uT§ET 
1 : 

2 ; 

It • 

4; 

5; 

6 ; 


11 
1 1 
7 
7 
7 
3 


1 

2 

3 

4 

5 

6 
7 
B 

10 

i| 

13 

If 

12 

9 

2 

4 

4 

2 

4 

4 

4 


5 

B 

5 


STOP 

e:no< of ekecution ^ 

CPU, riME!! 0.20 ELAPSED TIME; 


12 9 

U H 

12 11 


12 

12 


12 

12 


9 

9 

9 


12 11 
12 11 
3 

12 11 
12 11 


9 

9 


12 11 


13 11 

7 S 


13 

13 

13 

13 


9 

4 

3 

3 


13. 1 1 

ID 4 

3 5 

4 
3 

3 

4 
4 
7 


y 

5 

5 

5 

2 

3 


6 

2 

2 

3 

3 

3 

3 

1 

1 

1 

1 

1 

1 


5 

} 

1 

1 

1 

1 


1.18 


,038/BATCH^ 

CLGiTAJL AnotnaT 100 Is still iogg||^|D 4;o3'air 115100,15006411 
J,03. 31 user MANOD KUMAR IlSlOO , 1500641 _ 

Loggeid-Of I TTHIS at IfsSsslI ;on 13-^.Dtb:»B3: 

RUifclne:; 0;9o;oi, KCs;i6, eioooeof tise^s 0.;00;09- 



2 . 





1 


i:rC3M v»rsloni 102C2067) runnlEig Last saflgashcal 2445 in strean i 
pat' frsa DsfCC:LAST.CrUtl5103, 15^64] 
it'Oa.t tD: 3S!CC;LA5T.L0Gi:i5iaa, 1507643 
p ^aranatars 

Lae.:33j 71 ; 3 d Core;100p. Uniqae jYSB Paetart'j YES Qutput'UOuJG. 


uOGi:l 1 51DD/150064 /DEF£:R/;SP3aL:|yii/:n^S;: SO/CDRS: lOOP/LQCATSI "MiiNDJ 

33 1 7 I I r KANptJR 603 AC3) riYllS 
3tn:8:r jobs same pPM:24j 
3 23 2.9-;Sov-86 S'af 


30360 

folio 

lliu 

folio 

m 

ssip 

oEllo 


TY LA5T,f3R' 
3; 237 241 
0010 :i 
0020 :i 

0030 :i 

0040 :i 
'0050 :! 

■0060 

»007 0 Cl 
(0080 Ci 
>039 0 5i 
>0100 Ci 
>0110 C! 
)512D !!!i 

)0140 Ci 
0015 0 Cl 
00160 Ci 
00170 
00180 Ci 
00190 Ci 
00200 
00210 I 
00220 
00230 Cl 
00240 
00250 

00260 r 

niu 

00290 10 

00300. Ci 

0031 0 20 

00320 

00 330 
00340 


PROGRAM NAME:'*M30'TReS" 


.1 «# i. 


THIS PROGRAM MOOlFlJES tm FA’J.t^T TREE To TAKE CAPS- 3F 
XDR and JJOr GATES. XOR GArS: IS 5 H 3 WM BY ASSlGfJl'JG 13 
THE gate no. LiIARGEr THAN lOO.'NOT GATE IS INPiJ.Tl’EO^^BY 
assigning, it A! NgGATiyE NO.rHJStNOR and MAMD GATSP 
ARE ALSO TAKEN CARS- 3F.03.r?.7I GOMSISTS OF rRS| 

ING ONLY OR rANO GATES a-l3' BASIC' EVENTS EITHER NS-^ATEO 
OR UN-NEG.ArEO'.‘ , 

M OUT OF M kind OF GATES” AR'Si TAKEN CARE OF AT OOANl'in-- 
CAIION SrAG-E, 

THE program utilises DE, MpRiG-AN'S' LAW AND GATE: EOUlA^AueMC 


A|B=A-Bt A.B 


A-BsA.fB ^ A’t 3i“ Ai» 3: 


integer ck,gn(ioo#ioo),s,b,ad,or 

L«LEVEL OF TREE:. ORsNO.-^OF OR^^ GATES IN THE” TRESi 
MOsMAX, ORDER OF CUTSEIS («A)NrEO.' 
READC2§,1)L,0R,«^^ 

F0RMAT(lXr3ei3>lX)) 

INPUT FOR FAULir tRSE 
DO 10 I»l|100 

RE AD ( 2 2 f 2 5 , C GN ( I , J ) ,TPl , DM- 1 ) 

FORMATC 10014) , „ „ 

lFCGN(ia).E3.0)30rO 20 
BiBBfl 

continue 

B*=NO, OF BASIlC! EVENTS. 

s».o 

JPsO 

p’lao |sifig3 

DO 140 I si, 100 ^ ^ 
iFd.GT.aKOrO 130 
IF( J.GT# L3 GOTO 135 
lFtGN(I#J),Lr.400)GOT3 145; 

i?(5N(i;HS!:s:a.0)S.3r3 2S 
AOsAp+l 

AoSlEVEL pF.S'OB-BiASCH' ClONTA-LNlNG. THE: BASICI E.f SNT*” 

iFCAO.NE.LJGOro 29 

LSS'L 

^FtGJci,J,).€S.'JP)S-OTO 32* 

SSO 



I 


99. 
99 
99519 
99529 
99539 
9054 9 
90559 
905SD 
9 05 79. 

90589 

90590 
9.0590 
90510 
90629 
99539 
90649 
90659 

9958 9 

99719 

09720 

09730 

90740 

9878? 

09799 
99809 
90310 
99320. 
90830 
998 49 
00350 
90360 
903 7 9 
00380 
00390 
99900 
90910 
00920 

90939 

90940 

90959 

90960 
90979 
90989 
99999 
91000 
91010 

91029 

91030 
■ 1040 

1811 

0 I 8 I 8 

mu 

oil to 


32: 

35 


► mi W I*, .« " ’j , 


39 

32 


169 
17 9 


175 


49 


45 

48 


50 
7 0 


75 


145 

149 

lao 

''1 

5 

135 


iFcs.Sfs.o^soro 35 

CKsGf^Cl# J + l ) 

IP 1=1+1 
lP2=l+2 
JP1=J+1 
JP2=J+2 
LP1=L+1 
GN(I, J3=l 

IF(CK,NE,GMCljJPl))39r3 49 

lFCGMh,L.S).€5t0 3C^?P*6 

IFCGNCl,LS3..'iS:.0)L?P=6+l 

DO 80 j2sJ,LPP 

^^=LPP-J2 + J 

iFfM.LT. JP2)G9X3 82 

GNh>)=G^(I,Hil) 

CONTINUE 
GNtI,J + U=2 

Gi^CI. JP2)=.!l*3«Ci;, J.P2) 
8=B+l 

DO 170 I2SI.P2#8 

continue 

DO 175 J2=1#LP1 

CONTINUE . 

GN(IPI,JP,I3S4 

GN C IPl , J?2)=-:l ♦GN C i:Pl , JP29 
S“1 

GOTO 145 

I F ( GN ( I , LS ) . £9 . 0 J D P P s L 
ifCgn(i,ls5.#:o5dpP56+i 
DO 45 j2sJ>LPP 

IF(M^£t!jPZ)30T0 48 

GN(I.J+1)=2 

00^70 I2arp2.3! 

DO 60 j2sl,LH 
M=B««I2 + IP2 
Gn|m^J2|=GNCM-!1,J2) 

CgNTINUE: • 

DO 75 J2s,1,LP1 

’‘‘sGNClrJZ) 


lail^GMClPlf J?2i) 


costinoe 

GNapi/Ji , 

GNCIPl# JPi 

S*1 

continue 

CONTINUE 

RisTROciafllNG. OF tree; TO M ORDER *^ND LEVEL 0?. 
XREE,=' 

DO |60 1=1,6 

iFclJcifLiiy.'JE. 0330X0 OOO 

JlsL-Jtl 



3 


01125 

51130 

51140 

01150 

31 0 

01163 

300 

3117 3 


0110 0. 
OUOD 
0120 0 

oiilo 

01230 
012 43 
512'5 0 

230 

01260 

51270 

210 

01280 

200 

01290 

•«i| 

01305 

01310 

2!40 

01320 

mu 

250 

01350 

0i360 

260 

0137 0 



COfiTINOE 

FOLLOpNs rAt<;2s :;Ass 3 f 'lo-i, sor#»^i^o sates. 

DO 2DO 1=1, B 
•DO 2iO Js.i;l 

IFCGNCI, j;).GE.O)30rO 215 

ll i 5 ^ ^ t » J ) , 2 3 . 'J E . 0 3 3.3 13 2 30 

GNCI,J+l)s-il»3«(l , J + l) 

GOTO 210 
G«CI, J3=2 

GN { I - j+ p =-;i 3#< ( I , j:'»- 1 ) 

uOMTINUE 

CONTINUE 

OUTPUTTING r«.Sl MODIFIED FAJLiT TREE. 

WRITE(5,240)LitOR,M3 

FgRMAT(lX,3(i;3>lb) 

DO 260 1=1,8 

^811115 { 25 ol * LGM c I » J ) r J = l '» UfH ) 
FORMATClK,lO&r:4) 

CONTINUE 

STOP 

end 


,<r? f3R2 2-.'Da:t 
f3!23S27J 
1 115 

210 ! 

210 2 

0 

.Et bAST.fOR' 
l3:23S28i 

uoafing 

tUNKXcr uiAsr eixecution] 

2 1 10 
1 2 -1 

I 4 I 

I 2 2 

1 4 •12' 

STOP 

SINO* OF EKSCOTION , „ , 

CPU: Tl.MEr 5.04 ELAPSED TlMEJ 

sixi;r 


0,44 


.•KLJiOB/BArCH: 

tW-TAJL AnatWir iob is still 115100,1500643 3 

mo, 17 user manOJ kumar !i51 
U 03 ieji-bfl TflllS at' 3:23!3i- j- 
R|ntli»8;: . 0 :go:oO,,__KCS:iO, conoe.cf tiae.r 3:oO:ll 


ofsK'^Rftaa®: f34, writes 53 

3iA:rc3M vsrsIoD 102(2067) runolag LASr 5 8 
lOpat’ £ro» gSKCJL.A|T.CTL£15105.,r500||3 
O'oE'POit' t3l , DSKv J LAST » LOGi 1 51 05 , 1 5006 43 

fl5elo5:Sf':Oo* cortnoop uottasiyss 


an ao-'NJVrBS' 

50: 

alisfl 2 '4 47 in streai l- 


Rtstic’tiyss Output’S molog- 



i 


. LO 3.1 H 15130/150064 / D EF £R /S P 30 A tiiL/.T I vfES ’SO /] 
133 13 I I r fCANpOR 603AC3) 1^115 
CL'3.o;sp Ot.iair jobs same ?Pi4!2 4j 
03-25 sat' 

..ry F0R2 2'.<0AT 
t3:25M31 
3 3 5 

2 4 3 1 

2 4 S' 2 

2 4 S' 3 

2 201 3 4 

2 231 I 5 

2 2 31 1 6 

0 


:d R £.; 1 0 0 P / [i 0 3 a t s,:t i o /.?< e : “ m as 3 j k jm a = 


.€1 L AST. FOR' 

E3:25:i31 

Loadiog 

L»AST elxecutionl 


I, '*p«y • j 


2 

2 


2 

2 

2 

2 

2 

STOP 


4 

4 

4 

I 


3 

I' 

1 

2 

4 
2 
4 


0 

0 

0 

0 

0 

1 

2 

1 

2 


1 

1 

-4 

4 

5 

-5 

6 

••6 


S:fJ3' 3F EJCE'COTTON 

:iP-0, TIME Si 0.06 ELAPSED Tl HE S 0.36 
SiXI'T 

.CJOB/BATCE 

CLG.TAJfL Ahbtfteir.job is sfUl log3e3.-|a Jilalr £i5l00,i50064jj 
Ji03 13 J:S8:r HANOJ KUMAR l|5l 00 » 1530641 , 

Do ?' 3 e : d -: oC £ rtllis ,at 3 ; 25 S 15 , 3 ^ 23 - S 0 /;- Bj ; 

Root'liels 0sOO800» KCSs7r uonoBct tiiBed 3s‘03..O7 
D-lsC Reaisslei-f writes; 3 





A 38977 


TH- 

00^1 




Date sii|A38977 




is to be returned on the 


date last stamped. 




Mbtp- 1/187- (y) - Ku/yi ~fiur 


